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Executive  Summeury 


Many  military,  and  commercial,  communication  systems  operate  in  environ-* 
ments  containing  large  non-Gaussian  interference.  Tliis  interference  may  be 
hostile,  as  in  the  case  of  intentional  jamming,  or  friendly,  such  as  co-*channel 
interference.  Robust  locally  optimum  (LO)  detection  is  one  means  by  which 
interference  exhibiting  a  laxge  jammer-to-signal  ratio  {J/S)  can  be  mitigated. 
The  term  robust,  in  this  context,  indicates  that  prior  information  regarding 
the  channel  statistics  is  not  required.  Instead,  the  required  LO  nonlinear¬ 
ity  is  implemented  using  an  estimate  of  the  channel  noise  probability  density 
function  (pdf).  Thus,  the  robust  LO  detector  can  be  used  in  applications 
where  the  channel  interference  is  imknown  and  possibly  nonstationary. 

The  derivations  of  robust  LO  detection  algorithms  with  and  without 
memory,  and  their  performance  in  a  communications  system,  are  the  subject 
of  this  report.  There  aie  three  forms  of  LO  detection  considered:  (1)  detec¬ 
tion  with  memory,  (2)  detection  in  independent  and  identically  distributed 
(iid)  noise,  and  (3)  memoryless  detection.  A  number  of  algorithms  for  imple¬ 
menting  the  various  detector  structures  are  presented,  and  include:  (1)  the 
univariate  and  multivariate  histogram  algorithms,  (2)  the  M-interval  poly¬ 
nomial  approximation  (MIPA),  (3)  the  continuous  polynomial  approximation 
(CPA),  (4)  the  Fourier  series  approximation  (FSA),  and  (5)  the  univariate 
and  multivariate  kernel  algorithms.  In  all  of  these  methods,  observed  sam¬ 
ples  of  the  channel  noise  axe  used  to  estimate  the  noise  pdf,  from  which  an 
estimate  of  the  LO  nonlinearity  is  constructed.  However,  in  most  applica¬ 
tions  uncorrupted  noise  samples  are  imavailable  since  the  noise  is  corrupted 
by  the  information  signal.  A  common  solution  is  to  employ  the  large  JjS 
assumption  that  the  noise  pdf  is  approximately  equal  to  the  received  signal 
pdf.  However,  simulation  results  show  that  this  assumption  is  not  always 
valid,  even  in  high  J/S  environments. 
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To  determine  the  cause  of  the  observed  regions  of  poor  performance,  two 
“ideal”  detector  structures  were  examined:  (1)  the  globally  optimum  (GO), 
and  (2)  the  ideal  LO  (ILO)  detectors.  In  both  of  these  detector  algorithms, 
uncorrupted  noise  samples,  although  unavailable  in  practice,  were  provided 
and  used  to  estimate  the  required  noise  pdfs.  The  simulation  results  showed 
that  both  of  these  techniques  provided  “robust”  performance  over  a  wide 
range  of  jamming  scenarios.  It  can  be  inferred,  therefore,  that  robust  LO 
detection  techniques  will  be  extremely  useful  in  applications  characterized 
by  unknown,  nonstationary  interference,  provided  that  it  is  possible  to  ac¬ 
curately  estimate  the  noise  pdf.  As  a  result,  techniques  for  either  obtaining 
uncorrupted  estimates  of  the  noise  samples,  or  for  constructing  an  accurate 
estimate  of  the  noise  pdf  from  a  preliminary  estimate  of  the  received  sig¬ 
nal  pdf,  will  be  required  in  order  to  take  full  advantage  of  the  “robustness 
available  in  the  robust  LO  detector. 


8 


Chapter  1 
Introduction 


Many  military,  and  commercial,  communication  systems  operate  in  environ¬ 
ments  containing  large  non-Gaussian  interference.  This  interference  may  be 
hostile,  as  in  the  case  of  intentional  jamming,  or  friendly,  such  as  co-channel 
interference.  In  either  case,  a  performeince  degradation  will  likely  be  ob¬ 
served,  even  in  systems  employing  standard  spread  spectrum  (SS)  signaling 
techniques,  since  the  linear  receiver  structure  commonly  used  is  optimum 
only  in  the  ca^e  of  Gaussian  channel  noise. 

One  method  for  improving  system  performance  which  has  received  consid¬ 
erable  attention,  paxticuleirly  for  applications  in  high  jammer-to-signal  ratio 
{J/S)  environments,  is  locally  optimum  (LO)  detection  [l]-[4].  In  its  original 
form,  the  LO  detector  utilizes  a  nonlinearity  derived  cis  an  approximation  to 
the  optimum  mcasimum  likelihood  (ML)  detector.  In  many  scenarios,  the  LO 
detector  is  simpler  to  implement  than  the  corresponding  ML  detector,  and 
its  performance  asymptotically  approaches  that  of  the  ML  detector  as  the 
signal  becomes  small  relative  to  the  interference. 

One  drawback  to  the  original  form  of  the  LO  detector  is  that,  similar  to 
the  ML  detector,  it  requires  o  priori  knowledge  of  the  noise  probability  den¬ 
sity  function  (pdf)  to  implement  the  LO  nonlinearity.  In  most  applications 
of  interest,  however,  the  noise  statistics  are  unknown,  and  possibly  nonsta- 
tionaxy.  Thus,  recent  efforts  have  focused  on  robust  LO  detection  techniques 
[5]- [10],  in  which  the  required  LO  nonlinearity  is  constructed  using  estimates 
of  the  interference  statistics,  most  notably  the  pdf  and  its  derivative.  Since 
the  robust  LO  detector  implementation  is  based  on  estimation  techniques, 
it  can  adapt  to  the  unknown  noise,  and  modify  its  structure  as  the  noise 
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chcinges. 

The  derivations  of  robust  LO  detection  algorithms  with  and  without 
memory,  and  their  performance  in  a  communications  system,  are  the  sub¬ 
ject  of  this  report.  The  organization  of  the  report  is  as  follows.  The  deriva¬ 
tions  of  the  various  LO  detector  algorithms  for  quadrature  signaling  are  pre¬ 
sented  in  Chapter  2,  with  an  emphasis  on  direct  sequence  (DS)  SS  systems. 
Next,  a  number  of  memoryless  robust  LO  detector  methods  are  discussed 
in  Chapter  3.  These  methods  are  based  on  histogram,  polynomial,  Fourier 
series,  and  kernel  estimates  of  the  noise  magnitude  pdf.  Results  are  pre¬ 
sented  which  show  the  performance  characteristics  for  the  various  detector 
algorithms  when  subjected  to  a  single  continuous  wave  (CW)  jammer.  Algo¬ 
rithms  for  implementing  the  robust  LO  detector  with  memory  are  the  subject 
of  Chapter  4,  and  are  based  on  estimating  the  joint  multivariate  pdf  of  the 
channel  noise.  Results  are  also  presented  for  detection  in  various  CW  jam¬ 
mer  scenarios.  While  the  robust  LO  detectors  discussed  in  Chapters  3  and  4 
require  estimates  of  the  noise  pdf,  the  large  J /S  a.ssumption  that  the  noise 
pdf  is  approximately  equal  to  the  received  signal  pdf  is  used  in  simulation, 
since  uncorrupted  noise  observations  are  typically  unavailable  in  a  communi¬ 
cations  system.  The  effects  of  this  assumption  are  investigated  in  Chapter  5 
through  the  examination  of  two  “ideal”  detector  algorithms.  Finally,  the 
overall  conclusions  of  this  research  effort  eire  detailed  in  Chapter  6  of  this 
document. 
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Chapter  2 

Derivation  of  the  Robust  LO 
Detectors  for  Quadrature 
Signaling 


In  quadrature  signaling  techniques,  the  information  signal  is  transmitted  over 
two  orthogonal  channels,  namely  the  in-phase  (I)  and  quadrature  (Q)  chan¬ 
nels.  At  baseband,  these  signals  can  be  represented  as  a  complex  waveform 
given  by  [11] 

5m(t)  =  Si^{t)  +  jsQ^{t),  m  =  (2.1) 

where  and  are  T  second  duration  signals  transmitted  on  the 

I  and  Q  channels,  respectively.  M  denotes  the  number  of  different  possible 
signal  pairs.  Sampling  at  a  rate  /,  =  1/T,  produces  the  vector  s^,  written 
as 

Sm=S/„-hjSQ„  (2.2) 

where  s^,  s/„,  and  Sq^  axe  vectors  of  length  N,  eind  N  =  TfTa.  As  an 
example,  consider  basebmd  quadrature  phased  shift  keyed  (QPSK)  signaling. 
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The  possible  transmitted  signal  pairs  are 


l+jl 


SjM+jSQM  = 


1-;1 

-1+Jl 

-1-il 


,k  =  l,...,N 


(2.3) 


where  Eb  =  ^  Jo  \  Sm(t)  P  dt  is  the  energy  in  each  signal  period.  The 
corresponding  QPSK  signal  constellation  is  shown  in  Fig.  2.1. 


X  J¥- 

( 

X 

1 

1 

j  T 

X 

X 

Figure  2.1:  QPSK  signal  constellation 

The  derivation  of  the  LO  detectors  for  quadrature  signaling  is  as  follows: 
Let  the  received  signal  vector  of  length  N  be  given  by 

r/  +  jrq  =  s/„  +  n/  +  j{sQ^  +  nq)  (2.4) 

where  M  is  the  number  of  possible  transmitted  signal  pairs,  and  nj  and  nq 
are  the  I  and  Q  noise  vectors  (jamming  plus  backgroxmd  interference).  If 
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the  observed  value  of  the  received  signal  is  pj  +  jpq,  and  the  transmitted 
signal  and  noise  are  independent,  then  from  [5]  the  maximum  likelihood  (ML) 
detector  for  this  scenario  is  given  by: 


Choose  the  transmitted  signal  pair,  (s/„,Sq„),  which  maximizes 

ln-[/ii/n<j (Pj  —  ^ImiPq  ~~  SQm)]  , 


(2.5) 


where  /njn<j(^/j  Vq)  is  the  joint  probability  density  function  (pdf)  of  the  noise. 
Approximating  Eq.  (2.5)  using  first-order  Taylor  series  expansions  3rields  the 
following  result: 


lll[/n/n<j(P/  —  Sj„,  Pq  —  Sq„)]  «  lll[/n7n<j(P/,  Pq)] 

-  eL.  ln[/njn<j(P/,  Pq)] 

(2.6) 

Since  the  term  ln[/n/nQ(P/>PQ)]  is  constant  for  all  m,  the  ML  detector  can 
be  approximated  by: 


Choose  the  transmitted  signal  pair,  (s/„,Sq„),  which  maximizes 
KP)  =  KPl:Pq)  -  ~  T,k=i  S/„(fc)g^^ln[/n7nQ(Pj,PQ)] 

lll[/njiiQ(P7,  Pq)]]  (2-7) 


=  Cs/(P/,  Pq)  +  sl„gq{pj,  Pq)  , 


where  gi{Pi,  Pq)  =  [giM . . . giMf  >  Sq{Pi,Pq)  =  [9qM---9qM?'^ 
and 


9iM 

9qM 


Q 

^^lll[/n7nQ(P/,PQ)] 

Q 


e;^/njng(pj,pQ) 

fniHQiPi,  Pq) 
a^^fninqiPi,  Pq) 
fninqiPi,  Pq) 


(2.8) 
•  (2.9) 


The  functions  gi^{p)  and  gi^{p)  of  Eq.  (2.8)  and  Eq.  (2.9)  are  called  the  LO 
nonlinearities  with  memory,  and  Eq.  (2.7)  describes  the  LO  detector  with 
memory. 
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The  LO  detector  of  Eq.  (2.7)  can  be  greatly  simplified  if  the  noise  samples 
axe  independent  and  identically  distributed  (iid).  In  this  case,  the  joint  pdf 
of  the  channel  noise  becomes 

N 

—  U  fninQ{Vli^)>VQ{^))  •  (2.10) 

k=l 

The  LO  nonlinearities  of  Eq.  (2.8)  and  Eq.  (2.9)  then  become 
9ikip)  -  9i(pi,Pq)  =  —^HfnmQipiyPQ)]  = 

^Pl  fniriQ  {PliPQJ 

(2.11) 

^PQ  JntnqyPliPQ) 

(2.12) 

where  pi  =  pi(k)  and  pq  =  PQ{k)-  Thus,  the  2iV-variate  functions  of 
Eq.  (2.8)  and  Eq.  (2.9)  reduce  to  the  functions  of  2  variables  given  by 
Eq.  (2.11)  and  Eq.  (2.12).  Finally,  the  LO  detector  of  Eq.  (2.7)  becomes: 

Choose  the  treinsmitted  signal  pair,  (s/„,sq^),  which  maximizes 
KPIjPq)  =  -  [5u(A:)5^1n[/„^Q(p/(A:),pQ(A:))]  (2.13) 

+  '®Qm(^)apQ(jfe)  l^[/nriQ(/>/(^),  ^q(A:))]]  . 

Further  simplification  of  the  LO  detector  results  if  the  channel  noise  sam¬ 
ples  are  iid  and  if  the  joint  pdf  has  bivariate  radial  symmetry.  Under  this 
assumption,  the  joint  noise  pdf  cam  be  written  as 

/n(»7)/(27r77) ,  0  <  <  27r 

fnmqiVhVQ)  =  S  (2.14) 

0 ,  otherwise 

where  fn{v)  is  the  envelope  pdf  of  the  noise,  n  =  yjnj  +  nq  is  the  magnitude 
of  the  noise,  and  =  tan“^(n<j/n/)  is  the  phase  of  the  noise.  The  radial 
symmetry  assumption  is  valid  for  many  interference  sources  of  interest  as  they 
have  random  phase  angles.  Even  a  constant  frequency  waveform  will  have  a 
vector  that  rotates  at  a  uniform  rate  relative  to  the  transmitted  signal  vector. 
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and  is  therefore  equally  likely  at  any  angle  [10].  Under  this  assumption,  the 
LO  detector  of  Eq.  (2.13)  becomes  the  memoryless  LO  detector,  and  is  given 
by: 


Choose  the  transmitted  signal  pair,  (s/„,sq„),  which  maximizes 


KPi>Pq)  =  Er=i 

+5qm(%(/>(*))  sin(5r(fc))] 


(2.15) 


where 

g{p(k))  = 


d 


dp(k) 


ln[/„(p(fc))]  + 


p(k) 


imfMk))  1 

fn{p{k))  ’^pik)’ 


(2.16) 


p{k)  =  ^pjik)+  p^{k)  is  the  observed  magnitude  of  the  received  signal,  and 
6r(k)  =  tan“^(/ij<j(fc)/p/(A))  is  the  observed  phase. 

A  discussion  of  the  LO  nonlinearity  is  necesseiry  at  this  point.  In  deriving 
the  approximation  of  the  joint  noise  pdf  in  Eq.  (2.6),  an  importemt  assump¬ 
tion  was  made,  namely  the  small  signal  assumption.  Since  only  the  first-order 
Taylor  series  was  used  to  approximate  ln[/njn<,(-)]  Eq.  (2.6),  it  is  necessary 
that  the  higher  order  terms  of  the  Taylor  series  decay  to  zero.  These  higher 
order  terms  all  contain  powers  of  si^{k)  and  SQ^{k).  If  |  s/„(jfc)  l-C  1  and 
I  ^Qm{k)  |<  1,  then  [5/„(A:)]*’,  [sq„(^)F,  and  [5/„(A:)]*’45q„(A;)]'^  all  approach 
zero,  and  the  higher  order  terms  of  the  Taylor  series  become  negligible. 

The  small  signal  assumption  is  also  equivalent  to  the  large  jammer-to- 
signal  ratio  {JIS)  assumption,  i.e.,  that  the  noise  is  much  larger  than  the 
information  signal.  To  see  this,  assume  that  ]  n  1»|  s  |,  where  n  is  the 
random  noise  variable  and  s  is  the  information  signed.  Let  fi  =  nj A  emd 
s  =  sjA,  where  A  >  0  is  chosen  such  that  |  s  |C  1.  If  the  pdf  of  n  is  fnijl) 
and  the  pdf  of  n  is  fdfi),  then  fdfi)  =  A/„(An)  [12].  Note  that  if  t;  is  the 
observation  of  n,  then  rj  is  the  observation  of  n,  where  ^  =  jy/A.  The  LO 
nonlinearity  corresponding  to  fr,{-q)  is  g{p)  =  -f'MIUip)  +  1/^,  where  (') 
denotes  differentiation.  The  LO  nonlinearity  corresponding  to  fniji)  is  g{p)  = 
-fkip)lfn{p)  +  ^/P  =  -AfkW)!  fn{Ap)  -H  1//5  =  Ag{p),  where  p  =  pi  A. 
Since  g{p),  the  LO  nonlinearity  corresponding  to  the  large  J/S  assumption, 
and  g{p),  the  LO  nonlinearity  corresponding  to  the  small  signal  assumption, 
differ  by  only  a  constant,  the  resulting  LO  detector  structures  are  equivalent. 
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Thus,  the  Taylor  series  approximation  of  Eq.  (2.6)  is  also  valid  under  the  large 
Jj S  assumption  commonly  used  in  m^y  military  communications  scenarios. 

Of  particular  interest  in  this  report  is  the  LO  detection  of  QPSK  signals 
in  a  direct  sequence  (DS)  spread  spectrum  (SS)  communications  system.  In 
the  DSSS  system  shown  in  Fig.  2.2,  the  information  signal  in  the  interval, 
[iT,  (2+1)-T),  is  given  by 

d{t)  =  -  iT)  =  [dUi)  +  idQ„(0]n(t  -  iT)  (2.17) 

where  n(t)  is  a  unit  pulse  with  a  duration  of  T  seconds,  and  dm{i)  can  be 
one  of  four  (M  =  4)  possible  quadrature  signals: 

l+jl 
l-7l 

(2.18) 

-1+jl 
-l-jl 

The  information  signals  in  the  I  and  Q  channels  are  then  multiplied  by 
psuedo-noise  (PN)  “chipping  signals,”  c/(t)  and  CQ{t),  respectively,  as  given 

by 

=  /^  E  -  m  (219) 

V  ■‘c  fc 

where  Ec  =  2-^0"  I  c/(t)  +  jcqit)  dt  and  c/(fc)  =  ±1  and  cqijt)  =  ±1  each 
with  probability  1  /2.  The  resulting  transmitted  signal  during  the  interval 
in  the  QPSK  DSSS  system  is  ^ 

s{t)  =  si{t)  +  jsq{t)  =  X)  +  jdq„{i)cq{k)]  n(t  -  kTc) . 

V  ■‘■c  k=l 

(2.20) 

The  value  of  N,  the  length  of  the  signal  vectors,  is  given  hy  N  =  T fTc  and 
is  the  number  of  chips  per  information  symbol. 

^Note  that  dm{i)  is  constant  in  the  interval  [iT,  (i  +  1)  •  T). 
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Transmitter  Channel  Receiver 


Figure  2.2:  A  QPSK  DSSS  communications  system  with  LO  detection 


As  shown  in  Fig.  2.2,  the  sampled  version  of  the  received  signal,  r/(fc)  + 
jrq^k),  in  the  interval  [zT,  (z  +  l)r)  is  formed  by  sampling  the  received  signal 
at  Tc  seconds  yielding  the  following  expression 


ri{k)  +  jrQ{k)  =  J Y^u^i)ci{k)  +  ni{k)  +  j 


(2.21) 

where  njik)  +  jnQ{k)  is  the  (complex)  sampled  version  of  the  channel  noise. 
Let  the  observed  value  of  the  received  signal  be  Pi+jpq-  By  analogy  to  the 
derivations  of  the  LO  detectors  at  the  beginning  of  this  chapter,  and  removing 
any  constant  multipliers,  the  resulting  LO  detectors  for  the  QPSK  DSSS  sys¬ 
tem  are: 


Choose  the  transmitted  signed  pair,  {di^{i),dq^{i)),  which  maximizes 


N 


kpijpq)  —  ~ 


k=l 


HfnmQiPi,  Pq)] 


+  HfnjnciiPi,  Pq)]  (2.22) 


^One  can  edso  sample  at  higher  rates.  Here,  sampling  at  the  chip  rate  was  chosen  so  as 
to  simplify  the  notation.  The  simulation  results  in  subsequent  chapters,  however,  utilize 
multiple  samples  per  chip. 
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LO  detector  with  memory 
»  r  s 


KPi>Pq)  "  ^^m(*)c/(A:)^^^ln[/„^Q(/jj(fc),PQ(A:))] 

Q 

+  ^Qmi‘^)^Qi^)'Q^^^j^^^[fninQ{pl{k),  PQ{k))]  (2.23) 

LO  detector  for  iid  noise 

KPi^Pq)  =  S  {<^/m(0c/(%[p(*)]  cos[^,(fc)] 

A:=l 

-^dQ^{i)cQ{k)g[p{k)]  sin[flr(A;)]}  (2.24) 

memoryless  LO  detector 


The  LO  detectors  described  in  Eq.  (2.22)  to  Eq.  (2.24)  all  require  knowl¬ 
edge  of  the  noise  pdf.  However,  in  most  cases  of  interest  this  knowledge  is 
either  unavailable,  or  the  interference  is  nonstationary.  In  these  cases,  the 
LO  nonlinearities  must  be  estimated,  either  by  estimating  the  noise  pdf,  or 
by  estimating  the  nonlinearity  directly.  This  approach  of  robust  LO  detection 
is  the  subject  of  the  remainder  of  this  report. 
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Chapter  3 

Memoryless  LO  Detector 
Implementation  Techniques 


The  preceding  chapter  presented  the  derivations  of  a  number  of  LO  detectors. 
The  type  of  LO  detector  chosen  for  a  particular  application  is  primarily 
a  consequence  of  the  assumptions  made  concerning  the  noise  samples.  If 
the  noise  samples  cire  iid,  and  their  pdf  exhibits  radial  symmetry,  then  a 
memoryless  LO  detector,  such  as  that  of  Eq.  (2.15)  and  Eq.  (2.24)  may  be 
used.  Once  a  good  model  of  the  noise  is  obtained,  or  at  least  some  initiad 
assumptions  are  determined,  the  next  step  is  to  implement  the  corresponding 
LO  detector. 

In  the  military  environment,  since  most  intentional  jamming  techniques 
have  an  unknown  nature,  an  accurate  model  of  the  noise  pdf  is  usually  not 
known  a  •priori.  In  these  scenarios,  the  noise  pdf  must  be  estimated  from  the 
available  data  samples.  This  estimate  is  then  used  to  implement  an  approxi¬ 
mate  LO  nonlinearity.  If  the  uncorrupted  noise  samples  are  imavailable  (as  is 
the  usual  case  in  a  communications  system)  and  the  small-signal  assumption 
is  valid,  the  received  signal  scimples  may  be  used  in  many  cases  to  give  a 
crude  approximation  of  the  noise  pdf.  Thus,  various  univariate  pdf  estima¬ 
tion  techniques  and  their  use  in  implementing  the  LO  detector  are  discussed 
in  the  remainder  of  this  section. 
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3.1  Review  of  the  Histogram  Algorithm 

Using  the  histogram  as  a  means  of  implementing  the  memoryless  LO  non¬ 
linearity  has  been  extensively  studied  in  [5]  and  [6].  Recall  that  the  one¬ 
dimensional  memoryless  LO  nonlinearity  of  Eq.  (2.16)  can  be  written  as 

9{pi)  =  --^Hfnipi)]  +  -  •  (3.1) 

dpi  Pi 

To  approximate  g(pi)  using  the  histogram  algorithm,  the  first  step  is  to 
approximate  fniPi)-  Let  {7;^}  be  the  set  of  Q  observed  noise  magnitude 
samples,  and  divide  the  range  of  {77^}  into  K  intervals  or  “bins”,  Bk,  where 

Bk  =  {rii'  hk  <  rji  <  k  =  0, . . . ,  if  -  1}  ,  (3.2) 

The  values  {6fe},  with  bo  =  x^m  and  hK  =  ^max  >  are  called  the  “breakpoints” 
of  the  histogram  and  are  determined  heuristiceilly.  For  example,  the  break¬ 
points  may  be  chosen  so  that  the  bins  have  equal  width,  or  so  that  each  bin 
contains  the  same  number  of  observed  samples.  Next,  the  probability  of  a 
sample  being  in  each  bin,  P{Bk},  is  approximated  by  its  relative  frequency 

PW  =  A  S  IbM  (3.3) 

where  /^(y)  is  the  set  indicator  function,  i.e.,  /^(y)  =  1  for  y  6  -A  and 
0  otherwise.  Then,  the  histogram  estimate  of  fn{pi),  denoted  as  fn{pi),  is 
given  by 

/n(ft)=  E  (3.4) 

k=o  Pfc+l  - 

Next,  the  derivative  of  ln[/„(pi)]  must  be  computed.  If  it  is  assumed  that 
the  actual  pdf  fn{pi)  is  continuous,  then  the  impulses  that  arise  from  dif¬ 
ferentiating  Eq.  (3.4)  do  not  accurately  model  ^  lii[/n(/’i)]-  To  remedy  this 
dilemma  another  way  of  viewing  the  histogram  is  utilized.  One  can  think 
of  the  process  of  assigning  samples  to  a  bin  as  a  form  of  quantization.  In 
other  words,  all  pi  in  the  range  {bk  <  Pi  <  bk+i)  are  quantized  to  bk-  Then 
a  numerical  approximation  of  the  derivative  of  ln[/n(pt)]  evaluated  at  bk  is 
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used,  i.e., 


^  \pi=bH 


bk+1  —  £>*-1 

bk+i  -  fr*-i 


(3.5) 


M£>»42  -  ^Ai+i]  - 

£>fc+i  -  fcfc-1 


Note,  if  the  bins  have  equal  width  then  the  right  side  of  Eq.  (3.5)  reduces  to 
ln[P{Bfc4.i  }]  ^  jg  ^}jg  of  a  bin.  Finally,  the  value  of  the 

approximation  is  extended  from  bk  to  all  pi  in  the  range  {bk  <  pi  <  5fe+i), 
yielding  the  histogram  implementation  of  the  memoryless  LO  nonlinearity: 


^ — I 


iBkiPi)  +  —  • 

pi 


(3.6) 


Histogram  estimation  of  pdfs  usually  require  a  large  number  of  samples 
([14],  Ch.  3).  Thus  Q,  the  total  number  of  observed  samples  used  to  construct 
the  histogram,  may  be  much  larger  than  N,  the  number  of  samples  in  each 
signal  vector.  To  obteun  the  total  number  of  required  observed  samples,  the 
value  of  Q  should  be  chosen  such  that  Q  =  CN,  where  C  is  an  integer  chosen 
to  yield  enough  samples  to  construct  and  accurate  histogram.  The  C  vectors 
can  be  stored  and  all  the  available  samples  used  to  compute  the  histogram. 
Detection  can  then  be  performed  on  each  data  vector  using  the  resulting 
histogram  to  compute  the  approximate  nonlinearity  of  Eq.  (3.6). 


3.2  Review  of  the  MIPA  Algorithm 

The  M-interval  polynomial  approximation  (MIPA)  is  a  polynomieil-based 
method  for  estimating  the  noise  magnitude  pdf,  fn{v)  [5]  [9].  In  this  al¬ 
gorithm  the  domain  of  the  PDF  estimate  is  divided  into  M  intervals,  one 
for  each  approximating  polynomial.^  More  formally,  let  {rji}  be  the  set  of 
Q  observed  samples  of  the  noise  magnitude  random  variable,  n,  and  as  in 
the  histogram  method,  divide  the  range  of  into  M  intervals,  or  “bins,” 

^Note  that  in  this  context,  M,  the  number  of  intervals,  is  not  necessarily  the  same  as 
the  number  of  different  possible  transmitted  signal  pairs. 
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Bk  defined  by  the  bin  breakpoints  {5*}  as  in  Eq.  (3.2),  where  in  this  case 
M  =  K.  The  MIPA  estimate,  fniv)}  the  noise  pdf  can  be  written  as 


M-l  ^ 

fniv)  =  E  (3.7) 

i=0 

where 

A 

fiiv)  =  0,iQ  +  Oil  r)  +  ■  ‘  +  aipTf"  (3.8) 

is  the  approximating  pol3momial  for  the  bin,  p  is  the  MIPA  order,  and 
lA{y)  is  the  set  indicator  function.  The  MIPA  polynomial  coefficients  are 
obtained  by  first  constructing  a  polynomial  approximation  of  the  noise  pdf 
conditioned  on  the  event  that  the  observation  r}  lies  in  the  bin.  Recall 
that  the  conditional  noise,  faiv)}  for  the  bin  can  be  written  as 

fM  =  fnivl  hi  <V  <  hi+i)  =  p/L  N  (3.9) 

B[Oi  <7)  <  bi+i) 

where  P{bi  <r]<  bi+i)  is  the  probability  that  rj  is  in  the  interval  [6j,  6j+i). 
Let  the  approximating  polynomial,  faiv),  for  fciiv)  Ira-ve  the  form 

fciiv)  =  Cio  +  CiiV  +  '••  +  CipV^-  (3.10) 


The  coefficients  Cij  are  determined  by  minimizing  the  integrated  squared  error 
(ISE)  between  faiv)  fciiv)}  given  by 


ISE 


j  [fci{Tl)-L]^drj 

Jbi 


j 

hi 


6i+l 


l2 


faiv) Cij 

i=o 


dr} 


(3.11) 


Taking  the  derivative  of  the  ISE  with  respect  to  Cik,  setting  it  equal  to  zero, 
and  simplifying  yields  the  following  result: 


[  fMv^dr]  =1^7}^  Y^Cijn^dr},  A:  =  0, 

Jbi  Jbi 


,P- 


(3.12) 


The  term  on  the  left  side  of  Eq.  (3.12)  is  the  moment  of  faiv),  i.e.,  rrii^ 
[12].  Solving  the  integral  on  the  right  side  of  Eq.  (3.12)  and  writing  the  result 
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in  matrix  form  yields  the  solution  for  the  MIPA  coefficients  of  general  order: 


fci+l  —bj  ^1+1 
1  2 

2  3 

P+1 


(3.13) 


m,  =  Hj  X  Ci  (3.14) 

where  c  =  [cio  Cti  •  •  •  Cip]^  is  the  pol3momial  coefficient  vector,  m*  is  the 

vector  of  moments,  and  Hj  is  the  matrix  of  intervals.  From  this  formula  Cj 

can  be  computed  cis 

Ci  =  H,~^xmi.  (3.15) 

A  major  drawback  of  the  solution  for  Cj  in  Eq.  (3.15)  is  that  the  inversion 
of  Hi  is  computationally  intensive  and  must  be  performed  for  each  of  the  M 
intervals.  However,  a  significant  reduction  in  computation  can  be  achieved 
by  applying  the  linear  transform 

o»+i  —  Oi 

to  Eq.  (3.12).  This  expression  maps  the  interval  [6i,6i+i)  to  [0,1).  If  Xi  is 
the  random  variable  whose  outcome  is  a,  then  the  joint  pdf  of  Xi,  fxi{ot),  can 
be  written  as  fxi{oii)  =  (&i+i  —  h)  faiibi+i  —  bi)a  +  bi).  Thus,  the  system  of 
equations  in  Eq.  (3.12)  becomes 

=  «:  =  0 . p. 

Jo  Jo 

(3.17) 

Next,  examine  the  following  system  of  equations: 


/  fxi{ci)ct'‘da  =  I  dijcd'^'^da ,  A  =  0, . . .  ,p. 

j=0 


(3.18) 


The  solution  of  Eq.  (3.18)  for  d  =  [djo  dn  ...  djp]^  yields  the  MIPA  estimate, 
fxi{a),  for  fxiia),  namely 


fxi{a)  =  dio  +  diia  +  . . ,  +  dipO^ 


(3.19) 
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where 


di  =  X  fhi .  (3.20) 

The  vector  riii  conteiins  the  first  p  moments  of  Xi,  and  H  is  a  constant  matrix 
(for  all  i)  that  is  given  by 


1 


-  —  1 
2  *  •  •  p+l 


1 

2 


H  = 


1  _i_ 

3  • • •  p+2 


(3.21) 


-1- 

I  p+l  P+2 


-L- 
2p+l  J 


and  needs  to  be  inverted  only  once.  Comparing  Eq.  (3.17)  with  Eq.  (3.18),  it 
can  be  shown  that  the  coefficients  Cij  can  be  computed  from  a  linear  combina¬ 
tion  of  the  coefficients  dij  [].  If  L,-  is  the  corresponding  linear  transformation 
matrix,  then 

Ci  =  Li  X  di .  (3.22) 

Thus,  the  coefficients  aii  =  [uio  an  ...  ai^  for  the  overall  MIPA  pdf  estimate, 
/„(77),  in  Eq.  (3.7)  and  Eq.  (3.8)  can  be  found  as 

Hi  =  Li  X  X  ihi  •  P{hi  <T}  <  bi+i) .  (3.23) 


As  a  final  observation,  the  actual  vedues  for  the  first  p  moments  of  Xi  are 
usually  unavailable.  However,  they  can  be  estimated  from  the  available  data. 
Let  {pcj}  be  the  set  of  Qi  observations  of  the  noise  random  variable,  n,  that 
lie  in  the  bin.  If  mi^  is  the  A:*'*  moment  of  li,  then 


rhi^  = 


(3.24) 


can  be  used  to  estimate  mi,,.  Defining  riii  to  be  ihi  =  [1  rhij  . . .  mi^]^,  since 
mio  =  1,  then  the  MIPA  coefficients  in  Eq.  (3.23)  can  be  computed  as 


Once  the  MIPA  coefficient  vectors,  a^,  i  =  0, . . . ,  M— 1,  are  computed,  the 
resulting  MIPA  estimate,  ^(p),  of  the  memoryless  LO  nonlinearity  is  given 

by 


.f  .  ^  a-ii  +  2a, +  . . .  +  jaijp^  ^  +  . . .  +  pcippP  ^  1 

9{P)  =  -  2^  - ,  ,  ^  ■ - T— - hbiMiP)  +  T  • 

•  -  Clio  +  CLilP  + - +  CLipP^  p 


1=0 


3.3  The  CPA  Algorithm 

The  CPA  method  of  estimating  a  univariate  pdf  has  been  documented  in 
[7]  and  [8].  This  method,  like  the  MIPA  algorithm,  utilizes  a  concatena¬ 
tion  of  polynomials  to  construct  the  overall  pdf  estimate.  However,  the 
CPA  algorithm  also  provides  an  estimate  that  is  continuous  and  that  poss- 
eses  continuous  derivatives,  with  the  degree  of  continuity  determined  by  the 
approximation  order.  Four  types  of  CPA  algorithms  are  discussed  in  this  sec¬ 
tion:  the  basic  CPA,  the  CPA  with  linear  transform  (CPALT),  the  CPA  with 
Gaussian  tails  (CPAGT),  and  the  CPA  with  auxiliary  junction  (CPAUX). 

In  the  basic  CPA  algorithm,  an  estimate  of  the  cumulative  distribution 
function  (cdf)  corresponding  to  the  desired  pdf  is  formed.  This  estimate 
is  constructed  in  such  a  manner  so  as  to  have  continuous  first  and  second 
derivatives  at  all  bin  breakpoints.  The  result  is  a  pdf  estimate  which  has  a 
continuous  first  derivative,  a  property  that  is  desirable  when  implementing 
the  memoryless  LO  nonlinearity.  Formally,  given  {77^}  as  the  set  of  Q  ob¬ 
servations  of  the  noise  magnitude  random  variable,  divide  the  range  of  {77^} 
into  K  bins,  Bk,  with  breakpoints,  {6*},  defined  as  in  Eq.  (3.2).  Next,  a 
histogram  estimate  of  the  noise  cdf  is  constructed,  where  the  resulting  his¬ 
togram  values  at  the  breakpoints  60,  •  •  • ,  are  io, . . . ,  Yk,  with  10  =  0  and 
Yk  =  I-  Then,  the  basic  CPA  estimate,  Fii{rf),  of  the  noise  cdf,  Fn{v)j  is 
given  by 

=  E  yiMh,.k»)(v)  (3.2T) 

i=0 

where 

yi{T})  =  aio  +  anT]  + -■  +  aipif  (3.28) 

and  p  is  the  order  of  the  approximation.  The  value  of  p  is  determined  by 
the  number  of  continuity  constraints  required  in  the  algorithm.  For  the  basic 
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CPA  algoritlim  with  continuous  first  and  second  derivatives,  p  =  5,  and  is 
determined  by  using  the  constraints  that  Fn{r))  must  be  continuous  and  have 
continuous  first  and  second  derivatives  at  all  of  the  breakpoints  [7].  The 
coefficients  aij,  i  =  0,1,.. .  ,K  —  1  and  j  =  0, 1, . . . ,  p,  for  p  =  5  are  given  by 

[7] 
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where  Yi  is  the  value  of  the  histogram  cdf  at  the  breakpoint,  Y,'  is  the 
approximate  derivative  of  the  histogram  cdf  at  the  breakpoint,  and  Y"  is 
the  approximate  second  derivative  at  the  breakpoint.  The  values  for  Y/ 
and  Y"  are  computed  using  the  following  3-point  derivatives: 


>^^+1  - 


(3.30) 


and 


Y"  = 


Y!  —  V 

bi+i  —  bi-i 


(3.31) 


Note  that  in  Eq.  (3.30)  and  Eq.  (3.31),  the  expressions  for  Yq,  Y^,  Yq  and 
require  knowledge  of  b-2,  6-i,  ijc+z,  Y_2,  Y-i,  Yk+\,  and  Yk-+2. 

The  choice  of  these  additional  values  may  have  a  significant  impact  on  the 
resulting  CPA  estimate  [7].  Alternatively,  2-point  derivatives  may  be  used 
at  the  endpoints  to  eliminate  the  need  to  choose  additional  endpoint  values. 
In  any  event,  once  the  coefficients  {o^}  are  computed,  the  resulting  CPA 
estimate  of  the  memoryless  LO  nonlinearity  is 


2a;2  -f  Pcjap  -f  12ai4p^  -f  20at5p^  ^  i  1. 

^  Oil  +  2ai2P  +  SoiaP^  +  [6.A+i)  P  ^ 

The  CPALT  algorithm  is  similax  to  the  basic  CPA  method  except  that 
a  linear  transformation  is  performed  on  each  bin  of  the  histogram  cdf.  The 
strength  of  the  linear  transformation  is  that  it  greatly  reduces  the  complexity 
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of  the  equations  used  to  compute  the  CPA  coefficients.  In  particular,  the 
linear  transformation  is  given  by  the  pair  of  equations  [7]:  rji  =  r]  —  bi,  and 
Abi  =  bi+i  —  bi.  This  transformation  shifts  each  interval  to  the  origin.  The 
resulting  CPALT  estimate  of  the  noise  cdf  is 

K-l 

^n{v)  =  HVi)I[oAbi)(‘ni) ,  fji  =  r]-bi  (3.33) 

t=0 

where  the  polynomial  curve  for  each  interval  is 


Vi  (rji  )  =  Oio  +  Oil  rji  +  an  fj^  +  ••■  +  Ois  rjf  (3.34) 

and  {oLij}  are  the  transformed  polynomial  coefficients.  The  resulting  trans¬ 
formed  coefficients  for  p  =  5  are  computed  via 


■  1  Abi  Abf 

1  0  0 

0  1  2Abi 

0  1  0 

0  0  2 

0  0  2 


Abf  Abt  A5M  r  oio  1  r  Ti+i  • 

0  0  0  Oil  Vi 

3Abf  4Abf  5Abf  0,2  _ 

0  0  0  0,3  “  r/ 

6A6,  12A6?  20Abf  0,4 

0  0  0  J  L  J  [  r," 


(3.35) 


By  inspection,  0,0  =  T,,  a,i  =  Y-,  and  0,2  =  Y"  12.  Thus,  Eq.  (3.35)  can 
be  reduced  to  a  system  of  three  equations  in  three  unknowns,  resulting  in 
a  decrease  in  implementation  complexity.  Once  the  coefficients  {o,,}  axe 
computed,  the  CPALT  estimate  of  the  memoryless  LO  nonlinearity  is  given 
by 


N  20,2  d- 6a,3Pi  +  120,4/0,  +  20o,5P^  ^  1  _ _ _ 

M  =  -  2^  - — - r— T - -2^/  -3  ,  -4%.a6o(pO  +  -  >  (3.36) 

^  0,1  -(-  2o,2P,  +  3o,3pf  +  4o,4p?  -I-  5o,5p^  i  .  .a  /  \ 


where  pi  =  p  —  bi. 

The  CPALT  algorithm  may  be  further  simplified  by  invoking  the  fol¬ 
lowing  argument.  Recall  that  a  function  with  discontinuities  will  also  have 
discontinuous  derivatives.  Such  a  function  is  usually  an  unacceptable  esti¬ 
mate  for  most  continuous  pdfs  of  interest.  However,  the  CPALT  algorithm 
(as  well  as  the  other  CPA  algorithms)  utilizes  the  3-point  derivatives  given 
by  Eq.  (3.30)  and  Eq.  (3.31)  to  compute  the  estimated  deriavtives  at  the 
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breakpoints.  Thus,  an  auxiliary  function  can  be  used  in  place  of  the  initial 
cdf  estimate,  Fn{ri),  and  it  need  not  be  continuous  at  the  breakpoints.  How¬ 
ever,  the  auxiliary  function  must  satisfy  the  constraint  that  its  component 
first  and  second  derivatives  be  equal  at  the  breakpoints.  The  result  is  the 
CPA  with  auxiliary  function  method  of  estimating  the  noise  pdf  and  cor¬ 
responding  LO  nonlinearity.  Formally,  let  H{t])  be  the  axixiliary  function, 
where 

K-l 

Hiv)  =  S  yiiVi)hoAhi){Vi) )  Vi  =  T)-bi  (3.37) 

t=0 

and 

ViiVi)  =  d-  fii2rii  -}-  iSafjf  -I-  iSuTjf  .  (3.38) 

The  approximation  order  is  p  =  4  and  is  determined  by  the  constraints  on 
namely 


y'i{Vi)\fii=Abi 

=  y.-+i(^0h<=A&i  = 

yK%)lm=o 

=  yi-iivi)k=^o  = 

yi{vi)\vi=Abi 

=  yUv^)k=^b,  =  y:' 

p"(77i)hi=o 

=  yUvi)k=o  =  Yi, 

(3.39) 


The  coefficients  {Aj}j  i  =  0, . . . ,  if  -  1  and  j  =  1, . . . ,  4,  are  given  by 


1  2Abi  ZAbi  4A6i  ' 

■  Ai  ■ 

■  Y^  ' 

10  0  0 

YU 

0  2  6Abi  12Abi 

Y'' 

0  2  0  0 

(3.40) 
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Pi3  2A^ 

I3u  Abi(Y^’+Y^U)-2(Y;-YJ_,) 

L  4Ab? 

The  resulting  CPA  with  auxiliary  function  estimate  of  the  memoryless  LO 
nonlinearity  is 

->  /  N  +  G0i3pi  +  12/3i4Pi  T  f~  \  ^  ~  _  I 

(3.42) 

The  CPAGT  algorithm  is  particularly  useful  in  environments  where  Gaus¬ 
sian  noise  is  present.  Since  the  noise  observations  that  correspond  to  the  tails 
of  the  pdf  occur  infrequently,  it  is  difficult  to  estimate  these  regions  nonpara- 
metrically.  However,  if  a  Gaussian-like  tail  structure  is  assumed,  it  is  nec- 
esseiry  to  estimate  only  the  mean  and  variance,  since  these  two  parameters 
completely  characterize  the  Gaussian  pdf.  This  approach  is  the  foundation 
of  the  CPAGT  algorithm.  Beginning  with  a  histogram  approximation  of  the 
pdf  (not  the  cdf  as  in  the  basic  CPA  and  CPALT  methods),  the  CPAGT 
estimate,  /„(7/),  of  the  noise  pdf,  fn{v)  is 


jf-i 

fniv)  =  E  yi{v)I[bi.bi+i){v) 

i=0 

where,  in  this  case,  p  =  3  and 

=  aio  +  OiiT)  +  ai2Ti^  +  . 


(3.43) 


(3.44) 


The  CPA  coefficients  a,j, i=  1,2, ...  ,K  —  2  and  j  =  0, 1, . . .  ,p,  for  p  =  3  are 
given  by  [7] 


1  bi+i  bi+i  ii+i  1  r  Oio 

1  k  bj  Oil 

0  1  2bi+i  36|+i  ai2 

0  1  2bi  Zbf  J  L  o.i3 


(3.45) 
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However,  the  first  {bo  <  rj  <  bi)  and  last  (bK-i  <  V  <  ^k)  bins  of  the 
histogram  pdf  are  approximated  by  the  Gaussian  pdf 


where 
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(3.46) 


(3.47) 

(3.48) 


To  achieve  the  required  continuity  at  the  breeikpoints,  the  values  of  Yl, 
Yk-i,  and  are  computed  from  fg{r)).  The  resulting  CPAGT  estimate 
of  the  memoryless  LO  nonlinearity  is 


f  - + ^  h<p<bK-, 


M  =  { 


£rf* 

^2  7 


bo<  p  <  bi  OT 

bK-i  <  p  <bK 


(3.49) 

Note  that  a  linear  transformation  may  be  used  in  conjunction  with  this 
method  to  reduce  the  complexity  of  Eq.  (3.45). 


3.4  The  FSA  Algorithm 

This  section  presents  an  algorithm  for  estimating  pdfs  that  is  based  on  the 
Fourier  series  ([13],  Ch.  8).  Its  major  features  are  that:  (1)  the  approximation 
is  continuous  throughout  its  domain;  (2)  it  possesses  an  arbitrary  number  of 
continuous  derivatives;  and  (3)  no  quantization  process  is  required,  i.e.,  each 
data  sample  contributes  uniquely  to  the  overall  estimate. 

Given  the  set  of  Q  observed  noise  magnitude  samples  {77,}  with  mini¬ 
mum  value  rimin  and  maximum  value  rjmax,  the  Fourier  series  approximation 
(FSA)  of  the  noise  magnitude  pdf,  denoted  fn{ri),  is  given  by  the  following 
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expression: 

p 

^  {a*  cos(kujor])  +  bk  sm{ka}or})} 

for  77 min  ^  77  <  fjmax  (3.50) 

0  otherwise  . 

The  natural  frequency  of  the  FSA  is  given  by  wo  =  27r/r  where  T  = 
Vmax  ~  Vmin  Is  the  period.”  The  coefficients  {a^}  and  {6*}  are  determined 
by  minimizing  the  ISE  between  fn{v)  fn(v)y  actual  noise  magnitude 
pdf.  This  expression  is  written  as 

=  1:1:::  iUv)  -  MvWdrj 


Cfc  cos(A:a;o77) 


(3.51) 


p  _  12 

—  12  h  sin(fca;o77)  drj . 
fe=i 

DilFerentiating  Eq.  (3.51)  with  respect  to  oo,  a*,,  and  bk,  equating  each 
expression  to  zero,  and  noting  that  the  integration  of  the  sine  or  cosine 
functions  over  an  integer  multiple  of  its  period  is  equeil  to  zero  and  that 
fniv)  drj  =  1,*  yields  the  following  results: 


*It  is  assumed  that  the  support  length  of  the  actual  pdf  fn{T])  is  W*  -  ‘Hmin  ■  However, 
this  is  not  a  necessary  constraint,  since  truncation  can  be  used  to  approximate  a  pdf  with 
infinite  support. 
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Recalling  that  the  expectation  of  a  function  h{X)  of  the  random  variable  X 
is  given  by  E{h(X)}  =  h(x)  fx{x)dx,  Eq.  (3.53)  and  Eq.  (3.54)  can  be 
written  as 


“-fe  =  A:  =  i,...,p 

(3.55) 

Since  the  actual  pdf  of  n  is  unavailable  in  many  applications,  Eq.  (3.55) 
can  be  approximated  using  the  sample  averages,  resulting  in  the  estimated 
FSA  coefficients: 
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f]max  T/min  qf  j—j 


(3.56) 


Another  method  of  approximating  the  FSA  coefficients  is  discussed  in  [5]  and 
involves  first  estimating  fn{v)  ''^3-  ^  histogram  and  then  using  the  result  in 
Eq.  (3.53)  and  Eq.  (3.54).  If  the  value  of  the  histogram  at  the  breakpoint, 
bi,  is  Yi,  and  K  is  the  number  of  histogram  bins,  then  the  FSA  coefficients 
can  be  estimated  by 
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X)  COs{ku}o  k), 


^  max  V  min  i=o 


X  -  bi)sm{ku;obi), 


(3.57) 

A:  =  l,...,p 
k  =  l,...,p. 


To  implement  the  memoryless  LO  nonlinearity,  Eq.  (3.50)  and 
Eq.  (3.56)  are  used  directly  in  the  expression  for  g{pi)  given  in 
Eq.  (2.16).  The  resulting  estimate  of  g{pi)  is 


9{pi)  = 


X)  dfcA:a;o  sin(A;a;opj)  —  D  bkkuo  cos{ku}oPi) 
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P  P  /v  _ 
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Pi 


(3.58) 


32 


for  rjmin  <  pi  <  Vmax  and  0  otherwise. 

The  major  features  of  the  FSA  algorithm  can  be  seen  ihom 
Eq.  (3.50)  and  Eq.  (3.56).  First,  since  the  FSA  estimate  of  a  pdf  is  a  sum 
of  continuous  functions  the  overall  result  is  a  continuous  function.  Second, 
since  the  sine  and  cosine  functions  are  infinitely  differentiable,  the  result¬ 
ing  pdf  estimate  is  also.  Finally,  each  data  sample  contributes  uniquely  (no 
quantization)  to  the  FSA  estimate  through  Eq.  (3.56). 


3.5  The  Kernel  Algorithm 

Kernel  pdf  estimators  have  received  considerable  interest  in  the  past  few 
decades  ([14],  Ch.  6  and  references).  Like  the  histogram,  the  kernel  estima¬ 
tor  is  a  nonparametric  method,  i.e.,  it  does  not  assume  that  the  pdf  being 
estimated  has  a  particular  form.  Instead,  the  resulting  shape  of  the  pdf  esti¬ 
mate  is  solely  a  function  of  the  observed  data  samples.  The  kernel  estimator 
also  has  a  number  of  advantages  over  the  histogram:  (1)  it  is  continuous, 
(2)  it  may  possess  continuous  derivatives,  (3)  it  requires,  on  average,  less 
samples  to  form  an  estimate  than  does  the  histogram,  and  (4)  it  is  readily 
applicable  to  the  estimation  of  multivariate  pdfs  and  their  gradients.  Its 
major  disadvantage  is  that  in  many  cases  the  kernel  estimator  requires  con¬ 
siderably  more  computation  than  does  the  histogram.  However,  with  the 
availability  of  specialized  processors,  this  increased  computational  burden  is 
probably  not  an  issue. 

Since  the  subject  of  kernel  estimation  is  such  a  rich  area,  the  remainder 
of  this  section  provides  a  brief  overview  of  the  topic.  In  particular,  the 
kernel  estimators  for  univciriate  euad  multivariate  pdfs  will  be  presented,  with 
emphasis  placed  on  their  application  to  LO  detection. 

Let  {7?i}  be  the  set  of  Q  observed  samples  of  the  noise  magnitude  random 
variable  n.  The  kernel  estimator  of  fn{r]),  the  pdf  of  n,  may  be  written  as 

A(,)  =  ^P(^).  (3.59) 

The  kernel  window,  K{t),  determines  the  weighting  effect  that  the  samples 
{rji}  have  on  the  pdf  estimate  at  the  value  rj.  The  variable  h  is  known  as  the 
“smoothing  parameter”  and  determines  the  range  of  samples  {?;,}  that  will 
affect  the  pdf  estimate  at  the  value  rj. 
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Recall  that  the  memoryless  LO  nonlinearity  requires  the  first  derivative 
of  fn{v)-  Using  the  kernel  estimator,  this  can  be  approximated  as 


J, 


^  /r 


.  (3.60) 


As  can  be  seen,  the  kernel  estimator  can  provide  a  direct  approximation  of 
fn{v)j  which  relieves  the  need  to  use  numerical  differentiation  (as  in  the 
histogram  method). 

The  characteristics  of  the  kernel  pdf  estimator  of  Eq.  (3.59)  are  primarily 
determined  by  the  kernel  and  the  smoothing  parameter.  Kernels  usually 
satisfy  the  following  constraints: 

K{t)  >  0 

SK{t)dt  =  1  (3.61) 

StK{t)dt  =  0. 

Two  kernels  which  are  widely  used  are  (1)  the  Epanechnikov 

3 

K{t)  =  j(l  - 
and  (2)  the  Gaussian  kernel,  given  by 

K{t)  =  4=  ■  (3.63) 

v27r 

The  Epanechnikov  kernel  is  particularly  useful^since  it  minimizes  the  asymp¬ 
totic  mean  ISE  (AMISE)  between  /«(’?)  and  fn{ri).  Recalling  the  definition 
of  ISE,  e.g,  Eq.  (3.11)  or  Eq.  (3.51),  the  mean  ISE  (MISE)  is  given  by 

MISE  [/„(7,)]  =  E  {ISE  [/„(.))]}  ,  (3.64) 

and  the  AMISE  is  defined  to  be  the  main  terms  of  the  Taylor  series  approx¬ 
imation  of  the  MISE  ([14],  Ch.  2).  The  Gaussian  kernel,  on  the  other  hand, 
is  particularly  useful  in  situations  when  higher-order  derivatives  of  fniji)  are 
reqxiired.^ 

^For  cm  extensive  list  of  kernel  functions  see  reference  [14],  Chapter  6. 


kernel,  given 
(3.62) 
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The  smoothing  paxameter  h  also  alFects  the  AMISE  of  fniv)-  A  useful 
result  shown  in  ([14],  Ch.  6)  is  that  for  Gaussian  data  {77^}  and  a  Gaussian 
kernel,  the  value  of  h  which  minimizes  the  AMISE  is 

h  =  (3]  (3.65) 

where  is  the  variance  of  n.  Equation  (3.65)  can  be  a  useful  starting  point 
for  determining  a  value  of  h  suitable  for  the  given  application.  In  most  cases, 
the  value  of  h  will  need  to  be  decreased  as  the  data  becomes  “less  Gaussian.” 

The  simplest  way  to  apply  the  univariate  kernel  estimator  to  LO  detec¬ 
tion  is  by  direct  substitution  into  the  memoryless  LO  nonlinearity,  g{pi),  of 
Eq.  (2.16).  Thus,  the  kernel  estimate  of  g{pi)  is 


9{pi)  =  -77-7  +  -  = 

fn{pi)  Pi 


- — - 1 - , 


(3.66) 


More  sophisticated  forms  of  g{pi)  can  be  derived  if  characteristics  concerning 
the  noise  magnitude  pdf  axe  known  prior  to  detector  implementation. 


3.6  Simulation  Results  and  Comparison  of 
Techniques 

This  section  presents  the  simulation  results  for  the  memoryless  LO  detec¬ 
tor  implementations  discussed  in  Sections  3.1  through  3.5.  Results  are  pre¬ 
sented  for  LO  detection  in  both  a  standard  QPSK  system,  Eq.  (2.15)  and 
Eq.  (2.16),  and  a  QPSK  DSSS  system,  Eq.  (2.24).  In  implementing  the  LO 
nonlineaxities,  the  following  critical  assumption  is  made:  since  the  environ¬ 
ments  of  interest  exhibit  high  jammer-to-signal  {J/S)  ratios,  the  received 
signal  magnitude,  r,  is  approximately  equal  to  the  noise  magnitude,  n,  and 
thus  fr(-)  /„(•).  In  other  words,  the  received  magnitude  pdf  was  used  in 

place  of  the  noise  magnitude  pdf  in  the  implementation  of  the  memoryless 
LO  nonlineaxities.  Since  the  histogram,  MIPA,  and  histogram  FSA  methods 
were  examined  in  [5] ,  the  focus  of  this  section  is  primarily  on  the  various  CPA 
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algoritlinis  and  the  FSA  method  described  by  Eq.  (3.56).  Also,  analysis  of 
the  kernel  algorithm  is  postponed  until  Section  4.3. 

The  first  scenario  examined  was  the  probability  of  bit  error  (Pb)  perfor¬ 
mance  of  the  various  LO  nonlinearity  implementation  algorithms  in  a  stan¬ 
dard  QPSK  system  subjected  to  a  single  continuous  wave  (CW)  jammer.  The 
goal  of  analyzing  these  algorithms  in  a  QPSK  system  is  to  determine  parame¬ 
ter  settings  to  be  used  in  the  QPSK  DSSS  system.  Figure  3.1  through  Fig.  3.3 
illustrates  Pb  for  the  CPALT,  CPA  with  auxiliary  function,  and  CPAGT 
methods  for  various  values  of  fjfRa,  the  ratio  of  the  jammer’s  frequency  to 
the  information  signal’s  symbol  rate,  eind  N,  the  number  of  samples  per  sym¬ 
bol,  i.e.,  the  length  of  the  signal  vectors.'*  The  value  of  J/S  was  fixed  at  30 
dB  and  the  system  parameters  are  shown  in  Table  3.1.  As  can  be  seen  from 
the  graphs,  all  three  methods  exhibit  the  seime  performance  characteristics: 
(1)^6  decreases  as  fj/Rs  increases,  and  (2)  for  a  given  fjIRa,  Pb  remains 
relatively  constant  once  N  is  sufficiently  large,  approximately  N  >  20. 


System  Parameters  for  the  QPSK  System 

Pcirameter  Name 

Parameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Sampling  frequency  (in  Hz) 

1 

Samples  per  symbol  {N) 

5  to  50  by  5 

Number  of  CPA  bins  {K) 

10 

Number  of  samples  per  CPA  (Q) 

50,400 

Channel 

Parameters 

JjS  (in  dB) 

30 

filRs 

0.00496  to  0.496 
by  0.098208 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,008 

Table  3.1;  System  parameters  for  the  CPA  algorithms  in  a  QPSK  system: 
examination  of  Pb  relative  to  N  and  fjjR,. 

Next,  the  effects  of  K,  the  number  of  bins,  and  Q,  the  number  of  available 
samples  used  in  the  approximation  (i.e.,  the  number  of  received  magnitude 
samples),  on  Pb  was  examined  for  the  system  parameters  in  Table  3.2  and 

^The  results  shown  in  the  these  figures  were  compiled  from  only  a.  single  simulation  run. 
Thus,  these  graphs  should  be  used  for  the  purpose  of  determining  performance  trends,  and 
not  quantitative  values  of  Pj. 
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Figure  3.1:  Probability  of  bit  error  (1  CW  jammer  at  7/5=30  dB)  for  the 
LO  QPSK  simulation  -  CPA  with  linear  transform  implementation 


Figure  3.2:  Probability  of  bit  error  (1  CW  jammer  at  7/5=30  dB)  for  the 
LO  QPSK  simulation  -  CPA  with  auxiliary  function  implementation 

7 


3 


Figure  3.3;  Probability  of  bit  error  (1  CW  jeunmer  at  J/5=30  dB)  for  the 
LO  QPSK  simulation  -  CPA  with  Gaussian  tails  implementation 

JJS  =  30dB,  with  /j/iZa  =  0.496.®  The  Pi  plot  for  the  CPALT  algorithm, 
Fig.  3.4,  indicates  that  there  is  a  steady  decrease  in  Pi  as  both  K  and  Q 
are  increa.sed.  The  CPA  with  auxiliciry  function.  Fig.  3.5,  also  exhibits  a 
decrease  in  Pi  as  K  is  increased,  but  is  less  sensitive  to  varying  Q.  Finally, 
the  CPAGT  algorithm.  Fig.  3.6,  has  relatively  uniform  performance  for  all 
K  and  Q  exeunined  (provided  Q  is  sufficiently  large  for  the  given  K). 

Finally,  the  Pi  performance  of  the  various  CPA  algorithms  was  examined 
for  various  J/S  eind  K  values.  The  results  corresponding  to  /j/iZ*  =  0.496 
and  the  system  parauneters  in  Table  3.3  are  shown  in  Fig.  3.7  to  Fig.  3.9.® 
The  Pi  plots  indicate  that  all  the  methods  have  relatively  low  Pi  in  the 
region  2^dB  <J/S<  50dB,  provided  K  is  sufficiently  large,  and  also  in  the 
region  J/S  w  —lOdB.  However,  in  the  intermediate  region  of  —lOdB  < 
J/S  <  lOdB,  Pi  is  relatively  large.  While  the  reaison  for  this  increase  in 

®The  results  shown  in  the  these  figures  were  compiled  from  only  a  single  simidation  run. 
Thus,  these  graphs  should  be  used  for  the  purpose  of  determining  performance  trends,  and 
not  quantitative  values  of  Pj. 

®The  results  shown  in  the  these  figmes  were  compiled  from  an  average  of  three  simu¬ 
lation  runs. 


38 


System  Parameters  for  the  QPSK  System 

Pareuneter  Name 

Pewameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Sampling  frequency  (in  Hz) 

1 

Samples  per  symbol  (jV^) 

25 

Number  of  CPA  bins  (K) 

4  to  20  by  2 

Number  of  samples  per  CPA  (Q) 

5,000  to 
50,000  by  5,000 

Channel 

Parameters 

J/S  (in  dB) 

30 

hIRs 

0.496 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  cedculation 

10,000 

Table  3.2:  System  parameters  for  the  CPA  algorithms  in  a  QPSK  system: 
examination  of  Pb  relative  to  K  and  Q. 


Figure  3.4:  Probability  of  bit  error  (1  CW  jammer  at  J/5=30  dB)  for  the 
LO  QPSK  simulation  -  CPA  with  linear  transform  implementation 


Figure  3.5:  Probability  of  bit  error  (1  CW  jammer  at  J/S=Z0  dB)  for  the 
LO  QPSK  simulation  -  CPA  with  auxilieiry  function  implementation 


Figure  3.6:  Probability  of  bit  error  (1  CW  jammer  at  7/5=30  dB)  for  the 
LO  QPSK  simulation  -  CPA  with  Gaussian  tails  implementation 
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Pb  is  not  fully  known,  one  possible  explanation  is  that  under  the  large  J/S 
assumption  that  fr{-)  «  /n(‘))  LO  nonlinearity  attempts  to  mitigate  the 
dominant  component  of  the  interference.  For  20dB  <  J/S  <  50dB,  the  CW 
jammer  dominates,  while  for  J/S  —lOdB  the  Gaussian  background  noise 
is  dominant.  In  both  these  situations  the  LO  nonlinearity  is  able  to  mitigate 
the  dominant  interference,  and  the  effects  of  the  other  interferer  are  minima] 
However,  for  —lOdB  <  J/S  <  10d5,  neither  interferer  is  dominant,  so  the 
nonlinearity  must  attempt  to  mitigate  both,  with  a  resulting  increase  in  P5. 
A  further  investigation  must  be  performed  before  the  exact  cause  behind  this 
Pb  performance  characteristic  can  be  determined. 


System  Parameters  for  the  QPSK  System 

Parameter  Name 

Parameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Sampling  frequency  (in  Hz) 

1 

Samples  per  symbol  (AT) 

25 

Number  of  CPA  bins  {K) 

4  to  20  by  2 

Number  of  samples  per  CPA  (Q) 

5,000 

Channel 

Parameters 

J/S  (in  dB) 

-10  to  50  by  10 

fi/Rs 

0.496 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,000 

Table  3.3:  System  parameters  for  the  CPA  algorithms  in  a  QPSK  system: 
examination  of  Pb  relative  io  J/S  and  K. 

The  Pb  performance  of  the  FSA  implementation  of  the  memoryless  LO 
nonlinearity  in  a  QPSK  system  was  also  examined.  For  a  single  CW  jammer 
with  J/S  =  30dB,  the  Pb  relative  to  fj/R,  and  p,  the  FSA  order,  is  shown 
in  Fig.  3.10,  with  the  the  system  parameters  given  in  Table  3.4.^^  As  for 
the  CPA  algorithms,  Pb  decreases  for  the  FSA  algorithm  as  fj/Rg  increases. 
Also,  Pb  was  observed  to  be  lowest  in  the  region  5  <  p  <  20. 

With  the  information  garnered  in  the  analysis  of  the  memoryless  LO  non- 
linearities  in  a  QPSK  system,  a  Pb  performance  analysis  of  the  robust  LO 
nonlinearities  in  a  QPSK  DSSS  system  was  then  performed.  The  following 
nonlinearity  parameters  were  chosen  as  a  result  of  the  preceding  analysis. 

^The  results  shown  in  the  these  figures  were  compiled  from  an  average  of  three  simu¬ 
lation  runs. 
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Probability  of  Bit  Error 


Figure  3.7:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  simu 
lation  -  CPA  with  linear  transform  implementation 


J/S  in  dS 


Figure  3.8:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  simu 
lation  -  CPA  with  auxiliary  function  implementation 


System  Parameters  for  the  QPSK  System 

Parameter  Name 

Parameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Symbol  rate  {R3)  in  Hz 

0.05 

Samples  per  symbol  {N) 

20 

FSA  order  (p) 

2  to  20  by  2 

Number  of  samples  per  FSA  (Q) 

50,000 

Channel 

Parameters 

J/S  (in  dB) 

30 

filRs 

0.00496  to  0.496 
by  0.098208 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

5,000 

Table  3.4:  System  parameters  for  the  FSA  algorithm  in  a  QPSK  system: 
examination  of  Pb  relative  to  fj/Rt  and  p. 
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Figure  3.10:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK 
simulation  -  FSA  implementation 

For  the  CPA  algorithms:  K  =  10  and  Q  =  50, 400.  For  the  FSA  algorithm: 
p  =  10  and  20,  and  Q  =  50,400.  Using  a  processing  gain  =  TjTc  of  20 
(see  Chapter  2),  and  2  samples  per  chip  (resulting  in  A  =  2  x  20  =  40), 
the  performance  was  examined  for  the  various  robust  LO  nonlinearities  in 
a  QPSK  DSSS  system  subjected  to  a  single  CW  jammer  with  varying  JJS 
and  fjjRa,  with  the  system  parameters  provided  in  Table  3.5.®  The  sim¬ 
ulation  results  for  the  CPALT,  CPA  with  auxiliary  function,  and  CPAGT 
algorithms  are  shown  in  Fig.  3.11  to  Fig.  3.13.  Each  method  exhibited  the 
same  Pb  characteristic:  for  sufficiently  large  JjS  {J I S  >  lOdB)  Pb  decreases 
as  fj/Ra  increases.  Also,  the  CPAGT  algorithm  consistently  possessed  a 
Pb  lower  than  the  other  two  CPA  methods.  The  results  for  the  FSA  algo¬ 
rithm  with  p  =  10&20  are  shown  in  Fig.  3.14  and  Fig.  3.15.  This  algorithm 
also  exhibited  the  performance  trend  that  Pb  decreases  as  fj/Ra  increases. 
Finally,  similar  simulation  results  for  the  histogram,  histogram  FSA,  and 
second-order  MIPA  algorithms,  provided  from  [5],  are  shown  in  Fig.  3.16  to 

®The  results  shown  in  the  the  following  figures  were  compiled  from  an  average  of  three 
simulation  runs. 
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Fig.  3.18  for  compaxison  purposes. 


System  Parameters  for  the  QPSK  DSSS  System 

Parameter  Name 

Parameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Symbol  Rate  (22,)  in  Hz 

0.1 

Chip  rate  (i2c)  in  Hz 

2 

Samples  per  chip 

2 

Samples  per  symbol  (N) 

40 

Number  of  CPA  bins  (K) 

10 

FSA  Order  (p) 

10  and  20 

Number  of  samples  per 

CPA  and  FSA  (Q) 

50,400 

Channel 

Parameters 

J/S  (in  dB) 

0  to  50  by  10 

filRs 

igiraisii 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

5,040 

Table  3.5:  System  parameters  for  the  robust  memory  less  LO  nonlinearities 
in  a  QPSK  DSSS  system:  examination  of  Pi,  relative  to  JfS  amd  fj/Ra. 

From  an  examination  of  all  the  simulation  results,  a  number  of  obser¬ 
vations  can  be  made.  First,  while  some  LO  nonlineairity  implementation 
algorithms  may  exhibit  slightly  lower  Pi,  than  others,  in  general  they  all  have 
roughly  the  same  performance  characteristics.  Second,  and  probably  more 
importantly,  these  robust  techniques,  which  use  the  lairge  J/S  assumption 
that  /,(•)  w  /n(’)>  robust  as  would  be  desired.  In  particular,  for 

the  case  of  a  single  CW  jammer,  Pb  varies  greatly  as  fj/Ra  changes,  with  low 
Pi,  in  regions  of  high  fjIRa,  and  high  Pb  in  regions  of  low  fjIRa-  Thus,  the  LO 
detector  performs  poorly  when  the  frequency  of  the  CW  jammer  lies  near  the 
center  of  the  transmitted  signal’s  frequency  spectrum.  The  reason(s)  behind 
this  poor  performance,  and  a  discussion  of  possible  ways  to  improve  it,  are 
the  topics  of  Chapter  5. 
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Figure  3.11:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim¬ 
ulation  -  CPA  with  linear  transform  implementation 


Figure  3.12:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim¬ 
ulation  -  CPA  with  auxiliary  function  implementation 
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Figure  3.13:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim¬ 
ulation  -  CPA  with  Gaussian  tails  implementation 


Figure  3.14:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim¬ 
ulation  -  FSA  implementation  with  p  =  10 
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Figure  3.15:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim¬ 
ulation  -  FSA  implementation  with  p  =  20 


Figure  3.16:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim¬ 
ulation  -  histogram  implementation 
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j/SindB 


Figure  3.17:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim 
ulation  -  FSA/histogram  implementation 


Figure  3.18:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  DSSS  sim 
ulation  -  second-order  MIPA  implementation 


3.7  The  Direct  FSA  Implementation  of  the 
Memoryless  LO  Detector 


Given  a  few  assumptions,  a  Fourier  series  may  be  used  to  directly  estimate 
the  memoryless  LO  nonlinearity.  Ignoring  the  Ijp  term,  recall  that  the  mem¬ 
oryless  LO  nonlinearity  has  the  form  g{p)  =  —fn{p)  /  /n(p)-  To  derive  the 
FSA  estimate  of  g{p)  the  following  assumptions  are  made  concerning  the 
noise  pdf,  fniv)‘  fniv)  is  even  symmetric  about  the  origin  (and  thus  has  zero 
mean),  and  is  bounded  by  rimax  and  T^^tn,  where  rjmax  —  Begin  the 

derivation  by  approximating  g{p)  as 


M  = 


p  p 

^  +  E  cos(feu;o/9)  bk  sm{kujop) 

k=l  k=l 

for  T)  min  <P<  ^max 


0 


otherwise 


(3.67) 


where  T  =  rjrnax  —  Vmin  and  ujq  =  (27r)  /  T.  The  coefficients  {o^}  and 
{which  are  different  from  those  in  Eq.  (3.50))  are  chosen  so  as  to  minimize 
the  mean  squared  error  (MSE)  when  the  observation  consists  of  noise  only, 
i.e.,  r  =  n.  The  MSE  (denoted  by  the  function  J)  is  given  by: 

J  =  ^^{^(r)  -  5(0]^}  lr=n  =  7  i9{p)  -  9{p)?fn{p)dp 

— oo 


Vxnmx 

s 

Vmin 


9{P)  -  t 


p 

X;  ojt  cos{ku;op) 


X)  bk  sm{kujop) 

k=l 


fM  dp . 
(3.68) 


To  find  Co  which  minimizes  J,  the  expression  dJIdao  is  computed  and 

Vmnx 

equated  to  zero.  By  symmetry,  /  sm{ku;op)  fn{p)  dp  =  0  and  fn{Vmm)  = 

Vmin 

fn{v max),  resulting  in  the  following  equation: 

P 

^  £  “fe  E[cOs(A:a;on)]  =  fniVmin)  -  /n(77max)  =  0  .  (3.69) 

^  Jfc=l 

®Note:  These  cissumptions  are  valid  as  long  as  fniv)  bounded  and  is  symmetric  about 
its  mean  since  a  simple  linear  transformation  may  be  used  to  shift  the  pdf  to  the  origin. 
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Next,  the  expression  dJIdam  is  computed  for  m  7^  0  to  find  the  coeffi- 
cients  {a*}  which  minimize  J.  Noting  that  /  fn{p)  cos{rTuoop)dp  =  0  and 

**/"sin(ifca;op)cos(ma;op)/n(p)d/5  =  0,  the  resulting  system  of  equations  is: 

T?min 

p 

oq  E[c.o%{mu3on)]  2  a*,  {i;[cos((A:  -  m)a;on)] 
k=l 

+  £[cos((A:  +  m)a;on)]}  =  0,m  = 

(3.70) 

Finally,  J  is  minimized  with  respect  to  {6^}  by  computing  dJ jdhk  and  equat¬ 
ing  it  to  zero.  Using  the  facts  that  /n(7?inin)  =  /n(^max)  and  sin(Tna;o7/nim)  = 
sin(rm<;o?7m«)  =  0,^°  the  result  is  the  following  system  of  equations: 

i  f)  6ifc  {£?[cos((A:  -  m)a;on)]  —  i5[cos((fc -f- m)won)]} 
jfc=i 

=  /n(»?mm)sin(Tnu;o77imn)  “  /nCTymax)  sin(ma;o77nu«)  (3.71) 

-1-  rriu}QE[cos{rnu}Qn)] 


=  ma;oJ5[cos(77Won)] 


for  m  =  1, . . .  ,p.  To  simplify  notation,  define  the  vectors 


£L  —  [do  dj  . . .  dp] 

b  =  [h  h  ...  6p]^ 


E  {cos(u;on)} 

2  E  {cos(2u;on)} 
c  =  .  , 

{cos(pa;on)} 

'■°To  see  this,  recall  that  Wq  =  2irl{‘r]max  ~  »7mtn)  =  T^llmax  = 


(3.72) 


51 


eind  the  matrix 


1  —  E  {cos(2ti;on)} 


E  {cos(£>;on)} 

—E  {cos(3(i>o")} 


E  {cos[(P  -  l)a;on]} 
—E  {cos[(P  +  l)ti;on]} 


T  = 


E  {cOs(u»0Tl)} 

—E  {cos(3a<on)} 


1  —  E  {cos(4(i;on)} 


E  {cos[(P  -  2)u;on]} 
—E  {cos[(P  +  2)ii/on]} 


(3.73) 


E  {cos[(P  -  l)ta;on]} 
-E  {cos[(P  +  l)ti)on]} 


1  —  E  {cos(2Pti;on)} 


Then,  assuming  T  is  invertible  and  that  the  set  of  equations  described  by 
Eq.  (3.69)  and  Eq.  (3.70)  are  linearly  independent,  the  coefficients  {a^}  and 
{bk}  which  minimize  J  are  found  via  the  matrix  equations 


a  =  0 

b  =  2a;oT"’^c. 


(3.74) 


Thus  the  resulting  FSA  estimate  of  g{p)  is  given  by  (ignoring  the  1/p  term) 

M  sm{kuop)  (3.75) 

ifc=i 

for  J7inm  <  P  ^  Vmax  and  0  otherwise,  where  the  set  of  coefficients  {6*;}  is 
given  by  Eq.  (3.74). 

The  previous  derivations  for  the  FSA  estimate  of  a  pdf  (Section  3.4)  and 
of  the  memoryless  LO  nonlinearity  used  the  assumption  that  the  noise  pdf, 
fn{v)}  bounded.  However,  there  are  many  cases  of  interest  in  which  the 
support  of  fnijf)  is  unboimded,  e.g.,  the  Gaussian  pdf  [12],  Ch.  2).  In  these 
instances  the  FSA  algorithm  may  still  be  used  with  the  following  modifica¬ 
tions.  The  observed  received  signal  in  most  applications  will  have  a  maximum 
and  minimum  value.  If  fjmax  is  chosen  as  max(|  Pmax  |>  |  Pmin  |)>  then  the 
FSA  estimate  with  ljq  -  ir/ri^ax  will  encompass  all  values  of  the  received 
signal.  However,  the  actual  range  of  the  pdf,  (-oo,  oo),  can  be  used  to  com¬ 
pute  the  coefficients  via  Eq.  (3.55)  or  Eq.  (3.74).  While  the  result  will  not 
be  a  true  estimate  of  the  actual  unbounded  noise  pdf,  in  most  cases  it  will 
provide  a  suitable  approximation  given  the  fact  that  the  range  of  observed 
received  signal  values  is  boimded. 

The  FSA  algorithm  can  also  be  used  for  the  case  when  the  underlying 
pdf  is  unknown  and  knowledge  of  its  domain  (bounded  or  unbounded)  is 
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unavailable.  Provided  that  it  is  valid  to  assume  that  fn{v)  is  even-symmetric 
one  can  choose  rjmax  as  max(|  pzaax  |>  I  Pmin  !)•  The  required  coejfficients 
can  then  be  computed  using  Eq.  (3.56)  for  the  pdf  estimate,  or  Eq.  (3.74) 
(substituting  a  sample  average  for  the  expectation  operation)  for  a  direct 
estimate  of  the  LO  nonlinearity. 
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Chapter  4 

LO  Detector  with  Memory 
Implementation  Techniques 


The  implementation  of  a  robust  memoryless  LO  detector  was  the  subject  of 
Chapter  3.  In  this  chapter  two  of  these  techniques,  the  histogram  and  kernel 
algorithms,  are  extended  to  the  multidimensional  problem  of  implementing  a 
robust  LO  detector  with  memory.  As  before,  pdf  estimation  algorithms  are 
utilized  to  provide  the  freimework  by  which  the  required  LO  nonlinearities  of 
Eq.  (2.8)  and  Eq.  (2.9)  caji  be  estimated.  The  result  is  a  detector  that  cau 
be  used  in  an  unknown,  and  possible  changing,  interference  environment. 

As  will  be  seen,  the  algorithms  for  implementing  the  LO  nonlinearities 
with  memory  axe  more  complicated  than  their  memoryless  counterparts. 
However,  the  algorithms  with  memory  are  potentially  more  powerful  in  that 
they  utilize  the  interrelationship  between  data  samples,  whereas  the  mem¬ 
oryless  algorithms  process  each  Seimple  independently.  As  a  result,  the  LO 
detector  with  memory  cam  be  a  prudent  choice  in  environments  character¬ 
ized  by  highly  correlated  interference  samples,  such  as  channels  containing 
narrowband  jammers  or  noise  sources  modeled  as  auto  regressive  (AR)  pro¬ 
cesses. 


4.1  The  Histogram  Algorithm 

While  a  histogram  can  be  used  to  estimate  a  multivariate  pdf,  it  is  not 
paxticulaxly  suited  to  the  task  of  implementing  the  LO  nonlinearities  with 
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memory.  This  section  presents  the  multivariate  histogram  algorithm  and  the 
reasons  why  it  is  not  suited  to  the  task  of  implementing  the  LO  detector  with 
memory. 

Recall  that  the  I  and  Q  nonlineaxities  with  memory  of  Eq.  (2.8)  and 
Eq.  (2.9)  are  a  function  of  fniiiQ{Vi)VQ)i  the  multivariate  noise  pdf.  The 
histogram  estimate  of  a  multivariate  pdf  is  a  straightforward  (but  ardu¬ 
ous)  extension  of  the  histogram  algorithm  for  estimating  univariate  pdfs, 
discussed  in  Section  3.1.  Let  {rji.}  =  {[t;/^!)  -  ■  •  %(1V’)]^}  and  {rjQ.}  = 

c3'Ch  be  a  set  of  L  samples  of  the  I  and  Q  random 
noise  vectors,  n/  and  ng,  each  having  length  N.  Define  =  [nf  Hq]  as  the 


concatenated  random  noise  vector,  and  {rij}  = 

=  {[%(!)  •••  VliW  VQii^)  =  {[’?»(!)  •••  ’7t(21V)]}  as  the  set 

of  L  samples  of  n.  Thus,  /n/n<j(»7/>»7<j)  =  fniv)  is  a  21V-dimensional  pdf. 

Next,  divide  the  domain  of  the  dimension  of  fniv)  i^^to  K  regions 
using  the  breakpoints  bjk,  where  j  =  1, . . . , 2iV  and  k  =  0,. . .  ,K.  Thus,  the 
entire  2iV-dimensional  support  of  fniv)  is  divided  into  “hyper-bins,” 
Bki . defined  as 


=  {Vi  ■  bjkj  <  ViU)  <  bjikj+i)  for  j  =  1, . . . ,  2N 

and  kj  =  0, . . .  ,K  —  1}  . 


(4.1) 


Note  that  the  breakpoints  {bjk}  axe  chosen  heuristically,  e.g.,  such  that  each 
hyper-bin  ha^  the  same  width  in  each  dimension,  or  each  bin  contains  the 
same  number  of  sample  points. 

The  probability  of  a  sample  of  the  concatenated  noise  vector  lying  in  a 
given  bin,  P{Bki^...,k3tf}j  i®  approximated  by  its  relative  frequency. 


P{Bkx . »  P{Bki . kitf}  = 


1  ^ 

T  ■^^*1.— .*2W'  iVi) 

^  t=l 


(4.2) 


where  /^(y)  is  the  set  indicator  function.  Finally,  the  multivariate  histogram 
estimate  of  /n(*7)j  denoted  as  fniv)}  is  given  by 


K-l 


K-l 


fn(v)=E--  E 


1 


fel=0  k2jf=0  '^*1.— 


-P{Bk, . k2f,}Is^^ . k2jf(v)  (4.3) 
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where  the  normalizing  scale  factor,  is 

2N 

. •  (4-4) 

j=i 

At  this  point,  the  practicality  of  the  multivariate  histogram  estimator  of 
Eq.  (4.3)  should  be  discussed.  Suppose  that,  on  average,  it  is  desired  to 
have  P  samples  in  each  hyper-bin.  Then,  roughly  P  •  concatenated 
sample  vectors,  are  required.  The  result  is  a  required  total  number  of 
individual  noise  samples,  rii(j),  approximately  equal  to  2NPK^^  to  yield 
approximately  P  sample  points  in  each  hyper-bin.  For  example,  i{  K  = 
10  bins  per  dimension,  P  =  10  samples  per  hyper-bin,  and  N  =  2,  i.e., 
fnjnc}(Vi,VQ)  =  /n7(i)nr(2)n<j(i)nQ(2)(j7/(l),’7j(2),77Q(l),J7q(2)),  then  a  total  of 
approximately  400, 000  individual  noise  samples  are  required.  If  JV  =  10, 
however,  approximately  2  x  10^^  individual  noise  samples  are  required.  As  can 
be  seen,  a  linear  increase  in  the  length  of  the  I  and  Q  noise  vectors  results  in 
an  exponential  increase  in  the  number  of  individual  noise  data  samples.  Since 
in  most  situations  it  is  not  possible  to  obtain  the  required  large  number  of 
data  samples,  the  resulting  multivariate  histogram  estimate  will  be  poor.  An 
example  of  a  2-dimensional  histogreim  estimate  corresponding  to  iid  Gaussian 
I  and  Q  noise  data,  for  A"  =  6  and  L  =  10,000,  is  shown  in  Fig.  4.1. 

From  the  preceding  discussion,  it  is  evident  that  the  multivariate  his¬ 
togram  estimate  will  provide  an  extremely  crude  noise  pdf  estimate  at  best. 
Even  for  the  example  pdf  in  Fig.  4.1,  with  P  w  280  samples  per  hyper¬ 
bin,  the  resulting  estimate  is  poor.  In  addition,  the  numerical  differentiation 
techniques  required  to  produce  the  approximate  gradient  of  the  histogram 
estimate  will  introduce  further  error.  Thus,  a  multivariate  pdf  estimation 
technique  more  suited  to  the  implementation  of  the  LO  nonlinearities  with 
memory  is  presented  in  the  next  section. 


4.2  The  Kernel  Algorithm 


The  multivariate  kernel  estimator  of  the  noise  pdf  is  developed  a^  follows.  De- 

.r  _  r:T,Tj 


I  developed  a 

fine  the  concatenated  noise  vector  n^J  and  let 

=  {[77j(l)7/i(2)...77i(2Ar)]^},  i  =  be  a  set  of  L  samples  of 

the  random  noise  vector,  n.  Since  n  has  length  2N,  the  total  number  of 
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Figure  4.1:  Example  of  a  multivariate  histogram  estimate  of  a  pdf  -  iid 
Gaussian  I  and  Q  data  with  K  =  6  and  L  =  10,000. 

available  samples  is  =  2LN,  The  general  multivariate  kernel  estimator  is 
then  defined  as  ([14],  Ch.  6): 

Mv)  =  nVi 

The  window  iir(t)  is  now  a  multivariate  kernel  and  H  is  a  linecir  transform 
matrix. 

At  this  point  a  special  form  of  the  general  multiveiriate  kernel  estimator 
will  be  presented:  the  product  kernel  estimator.  The  product  kernel  estimate 
of  fn{v)  is  ([14],  Ch.  6): 

vU) 

(4.6) 

where  r]{j)  is  the  j*'*  argument  of  fn{v)  Vi{j)  is  the  sample  of  the 
observed  concatenated  noise  vector.  The  parameters  {hj}  are  the  smoothing 
parameters  corresponding  to  each  dimension  of  the  kernel  estimate.  The 
kernel  function  K{t)  is  the  same  univariate  kernel  discussed  in  Section  3.5. 


-  ViU) 
hj 


Mv)  =  /n(»7(l).  ",  Vii),  ",  ^(2A’))  = 


Lhi--- h2N  ^ 
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Recall  that  the  LO  detector  with  memory  requires  the  partial  derivatives 
of  the  multivariate  noise  pdf.  The  product  kernel  estimate  of  Q^fni'n)  is 

found  by  differentiating  fniv)  i^  Eq.  (4.6).  The  result  is 


Note  that  the  expression  in  Eq.  (4.7)  can  be  extended  to  partial  derivatives 
of  arbitrary  order. 

As  in  the  case  of  the  univariate  kernel  estimator,  the  product  kernel  es¬ 
timator  has  a  useful  expression  yielding  the  set  {/ij}  which  minimizes  the 
AMISE.  If  a  Gaussian  kernel  is  used  and  n  is  a  multivariate  Gaussian  vec¬ 
tor  with  independent  components,  then  the  value  of  hj  which  minimizes  the 
AMISE  of  fniv)  is  ([14],  Ch.  6) 


hi 


/  4  Nl/(2Ar+4) 

\2N  +  2) 


2^-1/(2JV+4) 


(4.8) 


where  <7j  is  the  variance  of  the  component  of  n.  Equation  (4.8)  provides 
a  starting  point  for  determining  {hj}  that  is  most  suitable  for  the  given 
application. 

The  product  kernel  estimate  can  be  used  to  approximate  the  LO  non¬ 
linearity  with  memory  by  substituting  the  kernel  estimates  of  fniv) 
e^fniv)  (9  =  Ij  •  ■  •  >  2i\r)  into  Eq.  (2.8)  and  Eq.  (2.9),  noting  that  fniv)  = 
/nrn<j(’7/(l)>  •  •  •  >  ViW,  ■■■,  9<3(-^))-  The  resulting  product  kernel  es¬ 

timate  of  gq{p)  =  —g^fn{p)/fn{p)  is: 
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The  techniques  employed  in  arriving  at  Eq.  (4.9)  can  be  used  to  develop  more 
sophisticated  forms  of  the  LO  detector  with  memory. 


4.3  Simulation  Results 


This  section  presents  the  simulation  results  for  the  product  kernel  implemen¬ 
tation  of  the  robust  LO  detector  with  memory.  Results  are  presented  for  LO 
detection,  with  and  without  memory,  in  a  standard  QPSK  system.  In  this 
case,  the  LO  detector  without  memory  refers  to  the  detector  of  Eq.  (2.13),  in 
which  successive  noise  samples  are  independent,  but  simultaneous  I  and  Q 
samples  are  correlated.  In  implementing  the  LO  nonlinearities,  the  following 
critical  assumption  is  made:  since  the  environments  of  interest  exhibit  high 
jammer-to-signal  {J j S)  ratios,  the  received  signal  vector,  vi-\-jvQ,  is  approx¬ 
imately  equal  to  the  noise  vector,  nj-fjnQ,  and  thus  /rjrQ(’)  fninqi')-  In 
other  words,  the  received  signal  pdf  was  used  in  place  of  the  noise  pdf  in  the 
implementation  of  the  LO  nonlinearities. 

The  first  scenario  examined  was  the  LO  detector  without  memory  in  a 
QPSK  system  subjected  to  a  single  CW  jammer.  The  P{,  performance  of 
the  product  kernel  algorithm  was  determined  relative  to  various  values  of 
J/S,  fjfRa,  and  B,  the  kernel  resolution.  The  parameter  B  is  related  to  the 
smoothing  parameters,  {hj},  by  the  following  expression: 


B  = 


Pjmax 


_ P Jmin 

hj 


(4.10) 


where  p,-  and  p,-  .  are  the  max:imum  and  minimum  observed  values  of 

r  jTTias  »  jTntn 

the  received  signal  in  the  dimension,  respectively.  The  resulting  Pb  plots 
are  shown  in  Fig.  4.2  to  Fig.  4.5,  with  the  system  parameters  provided  in 
Table  4.1.^  As  can  be  seen,  the  results  for  the  LO  detector  without  memory 
(but  with  correlated  I  and  Q  samples)  exhibits  a  similar  Pb  characteristic  as 
does  the  memoryless  LO  detector  algorithms,  whose  results  are  provided  in 
Section  3.6:  Pb  decreases  as  fjiRs  increases.  In  addition,  the  kernel  algorithm 
is  also  influenced  by  the  value  of  B.  For  good  performance  in  moderate  Jj  S 
regions  (lOdP  <  Jj S  <  20dB),  B  should  be  relatively  small  (10  <  B  <  20). 
However,  for  good  performance  in  large  J/S  regions  {20dB  <  J/S  <  AOdB), 

^The  results  shown  in  the  these  figures  were  compiled  from  an  average  of  five  simulation 
runs. 


59 


B  should  be  large  (20  <  B  <  40).  This  corresponds  to  the  discussion  in 
Section  3.6,  namely  that  hj  should  be  decreased  as  the  noise  environment 
becomes  “less  Gaussian.”  Thus,  as  J/S  increases,  the  interference  becomes 
“less  Gaussian,”  and  B  <x  1  /hj  must  be  increased. 


S3 

,rstem  Parameters  for  the  QPSK  System 

Parameter  Name 

Paxameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Symbol  rate  (P,)  in  Hz 

0.05 

Samples  per  symbol  {N) 

20 

Memory  length 

1 

Kernel  resolution  {B) 

10  to  40  by  10 

No.  of  samples  per  kernel  approx.  (Q) 

1,000 

Channel 

Parameters 

J/S  (in  dB) 

0  to  50  by  10 

filRs 

0.00496  to  0.496 
by  0.098208 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,000 

Table  4.1:  System  parameters  for  the  product  kernel  algorithm  (without 
memory)  in  a  QPSK  system:  examination  of  Pj,  relative  to  J/S,  fjIRs,  and 
B. 


The  next  scenario  examined  was  the  LO  detector  with  memory  in  a  QPSK 
system  subjected  to  a  single  CW  jammer.  Specifically,  the  value  of  N  chosen 
for  the  kernel  algorithm  was  iV  =  10  (see  Eq.  (4.9)).  The  plots  of  Pb  relative 
to  JJS,  fjIRs,  and  B  are  shown  in  Fig.  4.6  and  Fig.  4.7,  with  the  system 
parameters  provided  in  Table  4.2.^  As  before,  Pb  was  observed  to  decrease  as 
fjjRa  increased.  Also,  ss  J/S  increased,  an  increase  in  B  was  also  required 
to  maintain  a  good  vedue  of  Pb-  However,  the  observed  values  of  Pb  for 
the  LO  detector  with  memory  were  typically  larger  than  those  for  the  LO 
detector  without  memory.  The  most  probable  reason  for  this  phenomenon 
is  the  increase  in  error  associated  with  approximating  a  pdf  with  a  large 
dimensionality.  The  trade-off  between  memory  length  vs.  approximation 
complexity  will  need  to  be  addressed  in  future  research. 

From  the  analyses  made  in  this  section,  as  well  as  those  in  Section  3.6,  it 

^The  results  shown  in  the  these  figures  were  compiled  &om  an  average  of  five  simulation 
runs. 
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Figure  4.2:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  detector 
without  memory  simulation  -  B  =  10 


Figure  4.3:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  detector 
without  memory  simulation  -  B  =  20 


Figure  4.4:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  detector 
without  memory  simiilation  -  5  =  30 


Figure  4.5:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  detector 
without  memory  simulation  -  B  =  40 
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§3 

astern  Parameters  for  the  QPSK  System 

Parameter  Name 

Parameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Symbol  rate  (5,)  in  Hz 

0.05 

Samples  per  symbol  {N) 

10 

Memory  length 

10 

Kernel  resolution  (B) 

20  and  30 

No.  of  samples  per  kernel  approx.  (Q) 

10,000 

Channel 

Parameters 

J/S  (in  dB) 

0  to  50  by  10 

fi/Rs 

0.00496  to  0.496 
by  0.098208 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,000 

Table  4.2:  System  parameters  for  the  product  kernel  algorithm  (with  mem¬ 
ory)  in  a  QPSK  system:  examination  of  Pb  relative  to  JjS,  fjIRs,  and  B. 


Figure  4.6:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  detector 
with  memory  simulation  -  B  =  20 
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Figure  4.7:  Probability  of  bit  error  (1  CW  jammer)  for  the  LO  QPSK  detector 
with  memory  simulation  -  B  =  30 

can  be  concluded  that  the  LO  detectors  derived  using  the  large  J/S  assump¬ 
tion  that  the  received  signal  pdf  is  approximately  equal  to  the  noise  pdf  are 
not  as  robust  as  would  be  desired.  The  effect(s)  of  this  assumption,  as  well 
as  possible  methods  for  improving  performeince,  are  discussed  in  the  next 
chapter. 
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Chapter  5 

“Ideal”  Detection  Techniques 


The  simulation  aneilyses  of  the  LO  detector  techniques  discussed  in  Chap¬ 
ters  3  and  4  indicate  that  performance  varies  greatly  depending  on  the  jam¬ 
mer  characteristics.  In  particular,  for  the  case  of  a  single  CW  jammer,  all 
of  the  LO  detector  implementations  exhibited  the  same  trend  in  P^  in¬ 
creased  when  fj/Rs  decreased,  i.e.,  a  degradation  in  Pj,  was  observed  when 
the  CW  jammer  occurred  near  the  center  of  the  main  lobe  of  the  information 
signal’s  frequency  spectrum. 

However,  the  goal  of  the  robust  LO  detector  is  to  provide  reliable  de¬ 
tection  in  a  wide  range  of  possible  interference  scenarios.  In  an  effort  to 
determine  the  cause(s)  behind  the  observed  decrease  in  performance,  two 
ideal  detector  algorithms  were  examined:  the  globally  optimum  (GO)  and 
the  ideal  LO  (ILO)  detectors.  These  two  detector  schemes  are  ideal  in  the 
sense  that  uncorrupted  noise  samples  axe  assumed  to  be  available  at  the  re¬ 
ceiver.  Thus,  by  examining  these  two  detector  algorithms,  the  effects  of  the 
large  J/S  assumption  that  /n^ngC')  ^  /r/rgC')  can  be  further  investigated. 


5.1  The  Globally  Optimum  Detector 


The  GO  detector  for  a  standeird  QPSK  system  is  a  direct  implementation  of 
the  ML  detector  for  the  signal  and  additive  noise  scenario,  i.e.. 


Choose  the  transmitted  signal  pair,  (s/^,Sq„),  which  maximizes 
/njngCPj  —  S/„,Pq  —  Sq„)  , 


(5.1) 
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where  fmnqiVii  Vq)  is  the  joint  pdf  of  the  noise,  n/  +  jnQ,  and  Pi  +  jpq  is 
the  observed  value  of  the  received  signal  vector.  The  robust  GO  detector  can 
be  implemented  using  the  multivariate  kernel  estimator  discussed  in  Section 
4.2.  Define  n  as  the  concatenated  random  noise  vector,  =  nf  iiq  ,  and 

define  {rf[}  =  |  [77^  77qJ  }  to  be  the  set  of  L  vectors  formed  by  concatenating 
L  observed  samples  of  the  I  and  Q  channel  noise  vectors,  each  with  length 
N.  Then,  •  •  -  Vui^)  VQii^)  •  •  -  •  The  total  number  of 

available  noise  samples  is  thus  Q  =  2LN  Finally,  defining  the  following 
vectors,  =  pf  Pq  and  ,  and  using  the  result  in  Eq.  (4.6), 

the  product  kernel  implementation  of  the  robust  GO  detector  is 

Choose  Sto  which  maximizes 

(5.2) 

Up  -s^)  =  ZL  [n“.  K 

or,  more  explicitly. 

Choose  the  transmitted  signal  pair,(s/„,  Sq„),  which  maximizes 


In  Eq.  (5.2)  and  Eq.  (5.3),  hj  is  the  smoothing  parameter  in  the  sequence 
{hjj , . . . ,  hijf,  hq^ , . . . ,  hq^}  where  hi^  and  hq.  axe  the  smoothing  parameters 
for  the  dimension  of  the  I  and  Q  noise  vectors,  respectively;  r)i{j)  is  the 
sample  of  the  observed  concatenated  noise  vector,  77^,  and  correspondingly, 
rjliU)  VQiU)  samples  of  the  observed  I  and  Q  noise  vectors, 

respectively;  and  finally,  K{-)  is  any  one  of  a  number  of  kernel  windows 
discussed  in  Section  4.2  and  [14],  with  the  choice  of  kernel  determined  by  the 
requirements  specific  to  the  given  application. 

It  should  be  noted  that,  in  general,  the  GO  detector  cannot  be  imple¬ 
mented  in  practice  since  it  requires  uncorrupted  samples  of  the  channel  noise. 
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i.e.,  {riiiij),VQiU)}i  which  are  typically  unavailable  at  the  receiver  in  a  com¬ 
munications  system.  However,  the  analysis  of  this  detector  will  provide  useful 
information  concerning:  (1)  the  attainable  performajice  bounds  of  the  robust 
detectors,  (2)  the  effects  of  pdf  estimation  on  detector  implementation,  and 
(3)  the  implications,  and  inaccuracies,  associated  with  the  assumption  that 
^  friTqi')  ill  ^  high  J/S  environment. 


5.2  The  Ideal  LO  Detector 


Recall  from  Eq.  (2.7)  that  the  general  form  of  the  LO  detector  with  memory 
for  quadrature  signaling  is  given  by: 


Choose  the  possible  signal  petir,  (s/„,sq^),  which  maximizes: 
KPi^Pq)  =  [siM9lk{Pl,PQ)  +  5Qm(*)^c?k(p/,pQ)] 


where 


/  ^  d^f^i^giPi^PQ) 

QhiPuPq)  =  - (■  ^  . 

JnjnQiPi,  Pq) 


(5.5) 


and 


9Qk{PljPQ) 


dpq{k)f^I^QiPl’  Pq) 

fniHqiPi,  Pq) 


(5.6) 


are  the  ILO  nonlinearities  with  memory.  Using  the  result  given  in  Eq.  (4.9), 
and  the  definitions  in  the  previous  section,  the  product  kernel  implementa¬ 
tions  of  9i^{Pj,Pq)  and  gq^iPi^Pq)  are  given  by 


9Iu^PI<Pq) 


and 


Ef..| 

^  TTt  ( PiW-'^ZiW\  niv  /’p/O) 

1^  \ — ni  =  i 

1  i  #  *: 

EL  i 

[n”.  n,"..  IT  1 

''pQ(j)-n<5i(i)\ 

^  J 

} 

(5.7) 


9<3i,(PI’Pq) 
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where  (')  denotes  differentiation. 

As  with  the  GO  detector,  the  ILO  detector  is  not  usually  implementable 
since  it  requires  iincomipted  samples  of  the  channel  noise.  However,  the  ILO 
detector  can  be  utilized  to  determine  performance  bounds,  as  well  as  limi¬ 
tations  associated  with  LO  detector  implementations  that  use  the  received 
signal  samples  as  approximations  of  the  noise. 


5.3  Simulation  Results 

This  section  presents  the  simulation  results  for  the  GO  and  ILO  detectors 
discussed  in  Sections  5.1  and  5.2.  The  main  objective  behind  examining 
these  two  detector  structures  is  twofold:  (1)  to  determine  the  theoretical 
performance  limits  for  the  robust  LO  detector,  and  (2)  to  determine  the  effect 
on  Pb  of  the  large  J/S  assumption  that  the  noise  pdf  can  be  approximated  by 
the  received  signal  pdf. 

The  Pb  results  for  the  robust  GO  detector  in  a  standard  QPSK  system 
subjected  to  a  single  CW  jammer  are  shown  in  Fig.  5.1  to  Fig.  5.3,  with  the 
system  parameters  given  in  Table  5.1.^  In  these  figures,  Pb  is  shown  relative 
to  various  values  of  J/S,  fj/Rs,  and  B.  As  can  be  seen,  for  a  given  value 
of  J/S,  and  almost  any  fj/Rs,  Pb  Cein  be  reduced  to  useful  value  through  a 
judicious  choice  of  B.  In  particular,  for  low  J/S  {—20dB  <  J/S  <  lOdB),  a 
small  Veilue  of  H  (H  w  10)  is  required.  For  high  J/S  {20dB  <  J/S  <  4:0dB) 
a  larger  value  of  H  (20  <  H  <  30)  is  required.  These  results  support  the 
assertion  that  the  value  of  B  in  the  kernel  method  must  be  increased  as 
the  noise  becomes  “less  Gaussian.”  While  the  plots  for  the  GO  detector 
show  that  fj/Rs  has  little  effect  on  Pb  (for  appropriate  B  and  J/S  regions), 
the  same  was  not  observed  for  the  LO  detector  (see  Section  3.6).  Instead, 
Pb  for  the  LO  detector  increased  as  fj/Rs  decreased.  Thus,  the  observed 
performance  of  the  two  detectors  indicates  that  the  GO  detector  is  more 
robust  than  the  standard  LO  detector. 

Next,  the  performance  of  the  GO  detector  was  examined  for  the  case  of  a 

^The  results  shown  in  the  these  figures  were  compiled  from  only  a  single  simulation  run. 
Thus,  these  graphs  should  be  used  for  the  purpose  of  determining  performance  trends,  and 
not  qutintitative  values  of  Pj.  Also,  in  many  cases  the  resulting  calculated  value  of  Pj  was 
equal  to  zero.  However,  since  Pj  was  computed  using  only  1,000  sjrmbols,  a  value  of  zero 
actually  corresponds  to  a  P\,  on  the  order  of  10“®  or  less. 
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S: 

l^stem  Parameters  for  the  QPSK  System 

Parameter  Name 

Paxameter  Value 

Transmitter 

and 

Receiver 

Parameters 

Symbol  rate  (R*)  in  Hz 

0.05 

Samples  per  s3rmbol  [N) 

10 

Memory  length 

10 

Kernel  resolution  (R) 

10  to  30  by  10 

No.  of  samples  per  kernel  approx.  (Q) 

10,000 

Channel 

Parameters 

J/5  (in  dB) 

-20  to  50  by  10 

filRs 

0.00496  to  0.496 
by  0.098208 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,000 

Table  5.1:  System  parameters  for  the  GO  detector  in  a  QPSK  system  sub¬ 
jected  to  a  single  CW  jammer:  examination  of  Pb  relative  to  J/S,  fj/Rs,  and 
B. 


Figure  5.1:  Probability  of  bit  error  (1  CW  jammer)  for  the  GO  QPSK  de¬ 
tector  with  memory  simulation  -  S  =  10 


Figure  5.2:  Probability  of  bit  error  (1  CW  jammer)  for  the  GO  QPSK  de¬ 
tector  with  memory  simulation  -  B  =  20 


Figure  5.3:  Probability  of  bit  error  (1  CW  jammer)  for  the  GO  QPSK  de¬ 
tector  with  memory  simulation  -  B  =  30 
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QPSK  system  subjected  to  two  CW  jammers.  Specifically,  tbe  first  jammer 
was  fixed  at  J/5  =  30d5  and  fj/Rs  =  0.496,  with  the  second  jammer 
allowed  to  vary.  The  corresponding  Pj,  plots,  for  the  system  parameters  in 
Table  5.2,  are  provided  in  Fig.  5.4  to  Fig.  5.6.^  These  figures  show  that 
good  Pb  performance  can  be  attained  for  a  wide  variety  of  J/S  and  fj/R, 
second  jammer  configurations.  Previous  work  [5],  however,  has  shown  that 
the  standard  robust  LO  detector  does  not  perform  well  when  subjected  to 
two  CW  jammers,  except  when  one  has  a  low  J/S,  e.g.,  J/S  <  OdB. 


System  Parameters  for  the  QPSK  System 

Parameter  Neime 

Transmitter 

and 

Receiver 

Parameters 

S3rmbol  rate  (P,)  in  Hz 

0.05 

Samples  per  symbol  (N) 

10 

Memory  length 

10 

Kernel  resolution  (B) 

H|3Z£l:aEai 

No.  of  samples  per  kernel  approx.  (Q) 

10,000 

Channel 

Parameters 

J/S  (in  dB)  for  1*‘  jammer 

30 

fj/Rs  for  1**  jammer 

0.496 

J/S  (in  dB)  for  2”^  jammer 

-20  to  50  by  10 

fj/Rs  for  2"^  jammer 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,000 

Table  5.2:  System  parameters  for  the  GO  detector  in  a  QPSK  system  sub¬ 
jected  to  two  CW  jammers:  examination  of  Pb  relative  to  J/S,  fj/Rs,  and 
B. 


^The  results  shown  in  the  these  figures  were  compiled  from  only  a  single  simulation  run. 
Thus,  these  graphs  should  be  used  for  the  purpose  of  determining  performance  trends,  and 
not  quantitative  values  of  Pj.  Also,  in  many  cases  the  resulting  calculated  value  of  Pj  vras 
equed  to  zero.  However,  since  Pj  was  computed  using  only  1,000  symbols,  a  value  of  zero 
actually  corresponds  to  a  Pi  on  the  order  of  10~^  or  less. 
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Figure  5.4;  Probability  of  bit  error  (2  CW  jammers)  for  tbe  GO  QPSK 
detector  with  memory  simulation  -  B  =  10 


Figure  5.5:  Probability  of  bit  error  (2  CW  jammers)  for  the  GO  QPSK 
detector  with  memory  simulation  -  J9  =  20 
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Figure  5.6:  Probability  of  bit  error  (2  CW  jammers)  for  the  GO  QPSK 
detector  with  memory  simulation  -  5  =  30 

Finally,  the  performance  of  the  GO  detector  in  a  QPSK  system  subjected 
to  three  CW  jammers  weis  examined.  In  this  case,  the  first  jammer  was 
fixed  at  JfS  =  30dS  and  /y/iZ,  =  0.496,  the  second  at  J/5  =  ZddB  and 
fjiRa  =  0.496,  and  the  third  was  varied.  The  plots  of  Pj,  relative  to  JjS 
and  fjiRs  for  the  third  jammer,  as  well  as  various  values  of  B,  axe  provided 
in  Fig.  5.7  to  Fig.  5.9.^  The  corresponding  system  parameters  are  listed  in 
Table  5.3.  These  results,  as  before,  indicate  that  a  low  value  of  Pj,  can  be 
attained  for  a  wide  range  oi  JfS  and  fjfRs  settings  through  judicious  choice 
of  P. 

The  Pfe  performance  of  the  ideal  robust  LO  detector  was  also  examined. 
Specifically,  Pj,  for  the  ILO  detector  in  a  QPSK  system  subjected  to  a  single 
CW  jammer  is  plotted  in  Fig.  5.10  through  Fig.  5.12,  with  the  system  pa- 

^The  results  shown  in  the  these  figures  were  compiled  from  only  a  single  simulation  run. 
Thus,  these  graphs  should  be  used  for  the  purpose  of  determining  performcince  trends,  and 
not  quantitative  values  of  Ph-  Also,  in  many  cases  the  resulting  calculated  value  of  Pi  was 
equeJ  to  zero.  However,  since  Pt  was  computed  using  only  1,000  s3rmbols,  a  value  of  zero 
actually  corresponds  to  a  Pt  on  the  order  of  10“®  or  less. 
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System  Parameters  for  the  QPSK  System 


Parameter  Name  Parameter  Value 


Transmitter 

and 

Receiver 

Parameters 


Channel 

Parameters 


Symbol  rate  (iE,)  in  Hz 

0.05 

Samples  per  symbol  (N) 

10 

Memory  length 

10 

Kernel  resolution  (B) 

15  to  35  by  10 

No.  of  samples  per  kernel  approx.  {Q) 

10,000 

J/S  (in  dB)  for  1**  jammer 

30 

fj/Ra  for  1**  jammer 

0.496 

J/S  (in  dB)  for  2"^  jammer 

30 

fj/Rs  for  2"^  jammer 

0.406 

J/S  (in  dB)  for  3’’“^  jammer 

-20  to  50  by  10 

fj/Rs  for  S’""  jammer 

0.00456  to  0.456 
by  0.090288 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  ceilculation 

1,000 
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J3/S 

Fig\ire  5.9:  Probability  of  bit  error  (3  CW  jammers)  for  the  GO  QPSK 
detector  with  memory  simulation  -  5  =  35 

rameters  given  in  Table  5.4.'*  As  can  be  seen,  Pb  for  the  ideal  LO  detector  is 
similar  to  that  for  the  GO  detector  (see  Fig.  5.1  to  Fig.  5.3),  and  as  before,  an 
increase  in  JfS  requires  an  increase  in  B  in  order  to  yield  good  performance. 
However,  in  some  regions,  particularly  those  characterized  by  low  J/S,  the 
performance  of  the  ideal  LO  detector  is  poorer  than  that  of  the  GO  detector. 
This  is  most  likely  due  to  two  reasons:  (1)  because  of  computational  con¬ 
straints,  only  1,000  samples  were  used  in  each  ideal  LO  nonlinearity  estimate, 
whereas  10,000  samples  were  used  in  each  GO  detector  approximation,  and 
(2)  the  large  interfence  requirement  used  in  the  derivation  of  the  LO  detec¬ 
tor  is  not  met  when  J/S  is  small,  resulting  in  a  poor  approximation  to  the 
optimum,  or  GO,  detector  in  these  regions. 

The  resiilts  in  this  section  illustrate  a  number  of  important  points.  First, 
the  instances  of  large  Pb  for  the  LO  detector,  observed  in  Sections  3.6  and 

^The  results  shown  in  the  these  figures  were  compiled  from  only  a  single  simulation 
run.  Thus,  these  graphs  should  be  used  for  the  purpose  of  determining  performance 
trends,  and  not  quantitative  values  of  Pj.  Also,  in  many  cases  the  resulting  Pj  was  equal 
to  zero.  However,  since  Pj  was  computed  using  only  1,000  symbols,  a  value  of  zero  actueilly 
corresponds  to  a  Pj  on  the  order  of  10“^  or  less. 
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S3 

i^stem  Parameters  for  the  QPSK  System 

Parameter  Name 

Parameter  Value 

Transmitter 

and 

Receiver 

Parameters 

S3mibol  rate  (i2,)  in  Hz 

0.05 

Samples  per  symbol  {N) 

10 

Memory  length 

10 

Kernel  resolution  (B) 

10  to  30  by  10 

No.  of  samples  per  kernel  approx.  (Q) 

1,000 

Channel 

Parameters 

J/5  (in  dB) 

-20  to  50  by  10 

filRs 

0.00496  to  0.496 
by  0.098208 

Eb/No  (in  dB) 

10 

Symbols  used  for  Pb  calculation 

1,000 

Table  5.4:  System  parameters  for  the  ILO  detector  in  a  QPSK  system  sub¬ 
jected  to  a  single  CW  jammer:  examination  of  Pf,  relative  to  J/S,  fjlR$,  and 
B. 


Figure  5.10:  Probability  of  bit  error  (1  CW  jammer)  for  the  ideal  LO  QPSK 
detector  with  memory  simidation  -  B  =  10 
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Figure  5.11:  Probability  of  bit  error  (1  CW  jammer)  for  the  ideal  LO  QPSK 
detector  with  memory  simiilation  -  B  =  20 


Figure  5.12:  Probability  of  bit  error  (1  CW  jcimmer)  for  the  ideal  LO  QPSK 
detector  with  memory  simulation  -  5  =  30 
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4.3,  most  likely  can  not  be  attributed  to  error  in  the  nonlineaxity  approx¬ 
imations  since  these  same  estimation  techniques  were  used  in  the  GO  and 
ILO  detectors.  Second,  for  the  robust  LO  detector  in  a  system  subjected 
to  CW  jammers,  the  increase  in  Pb  corresponding  to  a  decrease  in  fj/Rs  is 
not  inherent  in  the  defining  expression,  Eq.  (2.7),  for  the  LO  detector,  as 
this  trend  weis  not  observed  in  the  Pb  results  for  the  ILO  detector.  Finally, 
the  large  J/S  assumption  that  the  noise  pdf  can  be  approximated  by  the  re¬ 
ceived  signal  pdf  is  the  most  likely  cause  of  the  LO  detector’s  regions  of  poor 
performance,  because  this  is  the  only  significant  difference  between  the  LO 
detector  and  the  GO/ILO  detector  algorithms.  Thus,  it  can  be  hypothesized 
that  improvements  in  the  performance  of  the  robust  LO  detector  cein  be  at¬ 
tained  by  constructing,  at  the  receiver,  good  estimates  of  the  noise  samples, 
or  an  accurate  estimate  of  the  noise  pdf. 
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Chapter  6 
Summary 


The  focus  of  this  report  has  been  robust  locally  optimum  (LO)  detection  in  a 
communications  system  subjected  to  interference  exhibiting  a  large  jammer- 
to-signal  ratio  {J/S).  The  term  robust,  in  this  context,  indicates  that  prior 
information  regarding  the  channel  statistics  is  not  required.  Instead,  the  LO 
nonlinearity  is  implemented  using  an  estimate  of  the  channel  noise  probability 
density  function  (pdf).  Thus,  the  robust  LO  detector  can  be  used  in  applica¬ 
tions  where  the  channel  interference  is  unknown  and  possibly  nonstationary. 

A  derivation  of  the  various  LO  detector  structures  was  provided  in  Chap¬ 
ter  2.  In  particular,  LO  detection  in  a  quadrature  signaling  environment 
was  examined,  with  emphasis  on  a  quadrature  phase  shift  keyed  (QPSK)  di¬ 
rect  sequence  (DS)  spread  spectrum  (SS)  system.  There  are  three  forms  of 
LO  detection:  (1)  detection  with  memory,  (2)  detection  in  independent  and 
identically  distributed  (iid)  noise,  and  (3)  memoryless  detection.  In  LO  de¬ 
tection  with  memory,  no  assumptions  concerning  the  noise  pdf  are  made, 
and  thus  the  entire  multivariate  pdf  is  used  in  constructing  the  LO  non¬ 
linearity.  LO  detection  in  iid  noise  utilizes  the  assumption  that  successive 
noise  seimples  are  iid,  but  with  simultaneous  in-phase  (I)  and  quadrature  (Q) 
samples  remaining  correlated.  This  detector  structure,  therefore,  utilizes  the 
two-dimensional  joint  pdf  of  the  I  and  Q  channel  noise  to  construct  the  LO 
nonlinearity.  Finally,  in  memoryless  detection  the  noise  is  assumed  to  be  iid, 
and  the  noise  pdf  is  assumed  to  possess  radial  symmetry.  In  this  case,  only 
the  univariate  pdf  of  the  noise  magnitude  is  required  to  implement  the  LO 
nonlinearity. 

With  the  derivations  of  the  various  LO  detector  structures  completed. 
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Chapter  3  focused  on  the  implementation  of  the  memoryless  LO  detector. 
Five  techniques,  which  utilize  estimates  of  the  noise  pdf  and  its  derivative, 
were  presented.  These  methods  are:  (1)  the  histogram  algorithm,  (2)  the 
M-interval  polynomial  approximation  (MIPA),  (3)  the  continuous  polynomial 
approximation  (CPA),  (4)  the  Fourier  series  approximation  (FSA),  and  (5) 
the  kernel  algorithm.  In  all  of  these  methods,  observed  samples  of  the  channel 
noise  are  used  to  estimate  the  noise  pdf  and  its  derivative,  from  which  an 
estimate  of  the  memoryless  LO  nonlinearity  is  constructed.  However,  in 
most  applications  uncorrupted  noise  samples  are  unavailable,  e.g,  the  noise 
is  corrupted  by  the  information  signal.  A  common  solution  is  to  employ  the 
large  Jj  S  assumption  that  the  noise  pdf  is  approximately  equal  to  the  received 
signal  pdf.  The  results  in  Section  3.6  illustrate  that  this  cissumption  is  not 
always  valid,  even  in  high  Jj S  environments. 

The  robust  LO  detector  with  memory  was  the  subject  of  Chapter  4.  Two 
methods  of  implementing  the  detector  were  presented:  (1)  the  histogram, 
and  (2)  the  kernel  algorithms.  The  histogram  algorithm  was  shown  to  be 
unsmtable  for  the  task  of  implementing  the  LO  nonlinearity  with  memory 
since  it  requires  an  extremely  large  number  of  observations  to  construct  a 
useful  pdf  estimate.  The  kernel  algorithm,  on  the  other  hand,  was  shown 
to  be  well-suited  to  this  task.  As  before,  since  uncorrupted  noise  samples 
are  usually  imavailable,  the  large  JjS  assumption  that  the  noise  pdf  can  be 
approximated  by  the  received  signal  pdf  was  utilized.  The  simulation  results 
in  Section  4.3  indicated  that,  as  for  the  memoryless  case,  this  assumption  is 
not  always  Vedid. 

The  simulation  results  for  the  robust  LO  detectors  with  and  without  mem¬ 
ory  implied  that  these  detectors  axe  not  as  robust  as  would  be  desired.  To 
determine  the  cause  of  the  observed  regions  of  poor  performance,  two  “ideal” 
detector  structures  were  examined  in  Chapter  5:  (1)  the  globally  optimum 
(GO),  and  (2)  the  ideal  LO  (ILO)  detectors.  In  both  of  these  detector  al¬ 
gorithms,  uncorrupted  noise  samples,  although  unavailable  in  practice,  were 
provided  and  used  to  estimate  the  required  noise  pdfs.  The  simulation  re¬ 
sults  showed  that  both  of  these  techniques  provided  “robust”  performance 
over  a  wide  range  of  jamming  scenarios.  Thus,  it  can  be  inferred  that  the 
poor  performance  regions  of  the  various  LO  detector  algorithms  were  caused 
by  the  use  of  the  large  J/S  assumption  that  the  noise  pdf  is  approximately 
equal  to  the  received  signal  pdf. 

In  conclusion,  robust  LO  detection  techniques  can  be  extremely  useful  in 
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applications  characterized  by  unknown,  nonstationary  interference,  provided 
that  an  accurate  estimate  of  the  noise  pdf  can  be  made.  Thus,  techniques 
for  either  obtaining  uncorrupted  estimates  of  the  noise  samples,  or  for  con¬ 
structing  an  accurate  estimate  of  the  noise  pdf  from  a  preliminary  estimate 
of  the  received  signal  pdf,  will  be  required  in  order  to  take  full  advantage  of 
the  “robustness”  available  in  the  robust  LO  detector. 
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MISSION 

OF 

ROME  LABORATORY 

Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Material 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence, 
reliability  science,  electro-magnetic  technology,  photonics,  signal 
processing,  and  computational  science. 

The  thrust  areas  of  technical  competence  include;  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


